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Summary 

y''  The  assumptions  and  criteria  used  in  existing  analyses  in  determining  the 
regions  and  transition  lines  of  pseudo-stationary  oblique-shock -wave  ref  lections  r-C 

have  been  re-examined  in  order  to  improve  the  agreement  between  experiments  and 
computed  data  for  regular  (RR) ,  single-Mach  (SMR) ,  complex-Mach  (CMR)  and  double- 
Mach  reflection  (DMR) . 

It  is  shown  that  the  relaxation  lengths  for  vibration  and  dissociation  deter¬ 
mine  whether  frozen  or  equilibrium  gas  transition  lines  are  applicable.  For  example, 
at  an  initial  temperature  of  300  K  and  a  pressure  of  15  torr  (where  much  previous 
work  was  done)  an  equilibrium-gas  analysis  would  not  be  required  for  shock  Mach 
numbers  Ms  <  9  in  N2 ,  Ms  <  6  in  02 ,  Ms  <  8  in  air  and  Ms  <  3  in  C02  . 

Yet,  the  available  experimental  data  in  N^,  C02  and  ■very  recent  results  for^3 
air^which”are  based  on  the  criterion  (consistent  with  relaxation  lengths)  of  the 
angle'ji*  between  the  incident  and  reflected  shock  wave,  do  not  conclusively  support 
the  frozen  or  equilibrium  gas  calculations  for  N"  and  air.  It  does  support  COj"  as 
an  equilibrium  gas  contrary  to  a  previous  conclusion  of  agreement  with  y  *  1.29. 

A  new  additional  and  necessary  criterion  for  the  transition  from  single  to 
complex  Mach  reflection  improves  the  agreement  between  analysis  and  experiment  and 
is  consistent  with  the  requirements  of  the  relaxation  length  and  the  angle  >  .  )<•  •  ■ 

However,  it  now  appears  that  a  more  accurate  criterion  is  required  for  the  boundary 
line  between  CMR  and  DMR. 

A  more  detailed  examination  of  the  boundary-layer-displq^ment  slope  at  the 
point  oi  regular  reflection  appears  to  eliminate  the  ^o-called  von  Neumann  paradox , 
and  explains  the  persistence  of  regular  reflection  below  the  transition  line  for 
the  occurrence  of  Mach  reflection. 

It  is  also  shown  that  at  the  triple  point  the  Mach  stem  can  vary  from  being 
perpendicular  to  the  wedge  surface  in  actual  experiments  by  as  much  as  -3.0°  to 
7.5°.  Consequently,  calculations  of  the  triple-point -trajectory  angle  y  on  the 
basis  that  the  stem  is  perpendicular  is  not  always  well  founded. 

It  is  verified  that  at  lower  shock  Mach  numbers  Ms  and  large  wedge  angles  9*, 
the  experimental  evidence  shows  that  the  transition  lines  for  SMR  *  CMR  and 
CMR  *■  DMR  converge  at  a  point  on  the  RR  *  MR  line,  contrary  to  a  previous  simplified 
analysis . 
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argon 

complex  Mach  reflection 
carbon  diox id© 

specific  heat  at  constant  pressure 
specific  heat  at  constant  volume 
double  Mach  reflection 
energy  of  vibrational  modes 
specific  enthalpy 
incident  shock  wave 


kink 

distance  from  wedge 
line 

distance  from  wedge 
horizontal  line 


comer  to  incident 


shock  wave  along  horizontal 


comer  to  kink  or 


second  triple  point  along 


Mach  number 


Mach  stems 
Mach  reflect*00 


incident  shock  wave  Mach  number 

now  Mach  number  in  region  1  relative  to 

flow  Mach  number  in  region  2  relative  to 
now  Mach  number  in  region  2  relative  to 
no  reflection 
nitrogen 


triple  point 
kink 

triple  point 


oxygen 

reflection  point 


pressure 
gas  constant 
reflected  shock  waves 


egular  reflection 
slipstreams 

single  Mach  reflection 

temperature 

triple  points 

velocity  of  incident  shock  wave 


/ 


1.  INTRODUCTION 

When  a  planar  shock  wave  collides  with  a  sharp 
compressive  corner,  four  different  types  of  reflec¬ 
tion  take  place  as  a  result  of  the  shock-wave 
reflection  and  the  deflection  of  the  shock-induced 
flow  'Refs.  1-4).  The  four  types  are  called  1  -eyular 
reflection  (RR),  single-Mach  reflection  (SMR), 
complex-Mach  reflection  (CMR)  and  double-Mach 
reflection  (DMR) .  The  latter  three  reflections  as 
a  group  are  called  Mach  reflection  (MR).  The  two 
types,  RR  and  SMR,  which  occur  also  in  steady 
supersonic  flows  were  first  noticed  by  E.  Mach  in 
1878  (Ref.  S) .  Much  later,  CMR  and  DMR  were  dis¬ 
covered  by  Smith  (Ref.  6)  in  1945  and  White  (Ref. 

7)  in  1951,  respectively,  while  working  in  shock- 
tube  flows.  The  four  types  are  illustrated  sche¬ 
matically  in  Fig.  1. 

The  simplest  type  of  reflection,  RR,  is  char¬ 
acterised  by  two  shock  waves,  an  incident  shock 
wave,  I,  and  a  reflected  shock  wave,  R,  and  their 
point  of  intersection,  P,  on  the  surface  of  the 
wedge.  Depending  on  the  wedge  angle  6W  and  the 
incident  shock  Mach  number  Ms,  the  intersection 
point  of  the  two  shock  waves,  P,  detaches  from  the 
wedge  surface  and  gives  rise  to  MR.  In  MR  a  new 
shock-wave  system  appears  including  Mach  stem  M, 
triple  point  T,  and  slipstream  (or  contact  surface) 
S.  In  SMR  the  curvature  of  the  reflected  shock 
wave  is  smooth.  In  CMR  a  kink  K  appears  in  the 
reflected  shock  wave  with  a  narrow  region  of  curva¬ 
ture  reversal.  In  order  for  the  kink  to  exist,  a 
band  of  compression  waves  must  converge  in  the 
region  behind  the  reflected  shock  wave  at  the  kink. 
This  band  is  usually  too  weak  to  be  seen  in  actual 
interferograms .  In  DMR  the  band  of  compression 
waves  converges  to  a  second  Mach  stem  M',  the  kink 
becomes  a  second  triple  point  T'  and  a  second  slip¬ 
stream  appears. 

Several  criteria  have  been  proposed  to  define 
the  transition  between  RR  and  MR.  von  Neumann 
(Ref.  8)  proposed  the  < detachment  criterion ,  where 
transition  takes  place  (see  Fig.  2)  when  the  wedge 
angle  4W  is  decreased  to  a  point  where  it  forces 
02  to  exceed  the  maximum  deflection  angle  9;>m 
(Ref.  9).  This  criterion  is  sometimes  called 
ih.  \  r  Neumann  criterion.  A  mechanical-equilibrium 
criterion  was  proposed  by  Henderson  and  Lozzi  (Ref. 
10)  on  the  basis  that  transition  can  only  occur 
smoothly  without  discontinuous  changes  in  pressure. 
This  criterion  can  be  formulated  as  01-02  =  ^3  *  0°. 
These  criteria  can  best  be  illustrated  by  using 
pressure-deflection  (p,  0)  shock  polars  as  shown  in 
Fig.  3.  The  reflected -shock  polar  R  corresponds  to 
transition  according  to  the  detachment  criterion. 

The  state  behind  the  reflected  shock  wave  jumps 
from  point  A,  which  corresponds  to  a  RR  state,  to 
point  B,  which  corresponds  to  a  MR  state.  Conse¬ 
quently,  there  is  a  pressure  jump  during  the  tran¬ 
sition.  The  polar  R*  corresponds  to  transition 
according  to  the  mechanical -equil ibrium  criterion. 
The  state  behind  the  reflected  shock  wave  comes 
along  the  p/po  axis  from  below  as  the  wedge  angle 
decreases  in  the  RR  region.  At  point  C  the  tran¬ 
sition  takes  place  and  RR  turns  smoothly  into  MR. 
After  the  transition,  the  state  behind  the  reflect¬ 
ed  shock  wave  goes  along  the  strong-shock  portion 
of  the  incident  shock  wave  I  to  the  right.  Further 
details  are  given  in  Ref.  1.  These  criteria  can  be 
expressed  in  another  way  as  follows.  The  detach¬ 
ment  criterion  is  the  limit  for  the  two-shock 


theory  to  have  solutions  and  the  mechanical- 
equilibrium  criterion  is  the  limit  for  the  three- 
shock  theory  to  have  solutions.  In  other  words, 
the  mechanical-equilibrium  criterion  corresponds 
to  x  =  0°  in  the  three-shock  theory  and  occurs  in 
stationary  supersonic  flows.  One  more  criterion 
was  suggested  by  Hornung  et  al  (Ref.  11).  Their 
criterion  is  called  the  sonic  criterion  because 
the  transition  takes  place  when  the  sonic  point  D 
(where  the  flow  behind  the  reflected  shock  wave  is 
sonic)  on  the  reflected  shock  wave  R"  coincides 
with  the  p/po  axis,  as  shown  in  Fig.  3.  This 
criterion,  however,  gives  a  transition  point  which 
is  so  close  to  the  point  of  detachment  that  it 
cannot  be  resolved  experimentally. 

The  first  precise  transit  ion -boundary  map  in 
a  (Ms-dw)  °r  (Ms~6w)  plane  based  on  the  two-shock 
theory,  the  three-shock  theory  and  the  transition 
criteria  described  above  was  given  by  Ben-Dor  and 
Glass  for  N2  and  Ar  (Refs.  1,  2).  Ando  and  Glass 
also  gave  a  transition-boundary  map  for  CO,  (Ref. 

13).  Ben-Dor  and  Glass  originally  suggested,  from 
a  comparison  between  experiments  and  calculations 
in  the  transition-boundary  maps,  that  in  the  range 
of  their  data  N2  should  be  treated  as  an  equili¬ 
brium  gas  and  Ar  as  a  frozen  gas  (>  =  1.667).  Ando 
and  Glass,  on  the  other  hand,  concluded  that  a 
fictitious  perfect -gas  model  with  ■,  =  1.290  for  C02 
is  the  most  appropriate  for  that  gas.  Although 
their  transition-boundary  maps  provide  very  useful 
information  for  practical  applications,  it  must  be 
concluded  that  any  apparent  agreement  with  experi¬ 
ment  is  accidental.  In  reality,  only  transition 
boundaries  based  on  a  frozen  flow  (>  =  1.667  for 
Ar  and  i  =  1.40  for  \2 ,  02  ,  air  and  C02 )  and  equi¬ 
librium  flow  are  consistent  with  relaxation  length 
concepts.  The  actual  experimental  data  docs  not 
agree  with  either  model  in  transition-boundary  maps. 
Consequently,  it  must  be  concluded  that  better 
criteria  for  the  SMR  *  CMR  and  CMR  *  DMR  transitions 
must  be  found,  which  would  hold  for  all  gases  tested 
to  date. 

A  much  better  and  fundamental  assessment  can 
now  be  given,  whether  a  flow  is  in  equilibrium  or 
not,  through  the  appropriate  relaxation  length  as 
represented  by  the  angle  i  between  the  incident  and 
reflected  waves  rather  than  by  the  secondary  stan¬ 
dard  of  agreement  with  Ms-Pw  plots.  A  good  case  in 
point  is  the  conclusion  by  Ando  and  Glass  (Ref.  13) 
that  the  fictitious  perfect-gas  (>  =  1.290)  model 
for  the  Ms-(*w  plot  is  the  one  that  agrees  best  for 
C02  over  the  range  1  <  Ms  <  10.  However,  a  detailed 
examination  of  6  over  this  range  shows  that  the  flow 
is  in  vibrational  equilibrium  because  the  vibration¬ 
al  relaxation  lengths  for  Ms  3  are  all  too  small 
(<  1  mm).  On  the  other  hand,  for  Ar,  N2 ,  02  and 
air,  this  shock-tube  data,  by  and  large,  can  be 
considered  as  consistent  with  perfect-gas  (frozen) 
states.  The  above  leaves  little  doubt  that  even 
more  precise  SMR  *  CMR  and  CMR  *  DMR  criteria  are 
still  required.  A  new  necessary  (but  not  sufficient) 
condition  for  SMR  *  CMR  transition  is  presented 
which  improves  the  agreement  with  experiment  and  is 
consistent  with  relaxation  lengths.  It  is  also 
shown  that  the  experimental  data  (Refs.  13-lb)  do 
not  agree  with  the  CMR  *  DMR  transition  line  at 
lower  Ms  and  large  Pw  In  fact,  the  two  lines  merge 
where  the  SMR  *  CMR  transition  line  cu*s  the  RR  *  MR 
line.  Basically,  this  arises  from  the  fact  that  the 
distance  between  the  two  triple  points  T  and  T'  does 
not  remain  constant  but  vanishes  at  the  point  of 
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intersection  of  the  foregoing  three  transition 
lines. 

Fot  nonstationaTy  shock  reflections,  experi¬ 
mental  results  show  that  transition  occurs  in 
accordance  with  the  detachment  criterion.  Experi¬ 
ments  also  show  that  RR  persists  beyond  the  limit 
of  'he  detachment  criterion.  This  persistence  has 
been  called  the  von  Neumann  paradox.  Homung  et 
al  (Ref.  11)  suggest  that  it  may  be  explained  in 
terms  of  the  displacement  thickness  of  the  boundary 
layer  on  the  wedge  surface.  Ben-Dor  et  al  (Ref.  12) 
showed  that  a  persistent  hysteresis  loop  exists  in 
the  RR  *  MR  transition  in  experiments  which  depend 
on  whether  6W  is  increased  or  decreased  for  a  given 
Ms.  The  reason  for  this  persistence  is  still  un¬ 
clear  and  will  be  discussed  in  Section  3.4.  Once 
RR  terminates,  three  different  types  of  Mach  re¬ 
flection  can  occur  in  nonstationary  flows.  White 
(Ref.  7)  noticed  that  CMR  and  DMR  will  occur  when 
the  flow  behind  the  reflected  shock  wave  becomes 
supersonic  in  a  frc.me  of  reference  attached  to  the 
triple  point  T  (i.e.,  M2T  '  1,  Fig.  2).  Henderson 
and  Lozzi  (Ref.  10)  suggest  that  a  band  of  compres¬ 
sion  waves  must  exist  in  a  CMR  and  these  compression 
waves  converge  to  a  shock  wave,  the  second  Mach 
stem,  to  form  a  DMR  when  M2k  >  1.  The  criteria  for 
transition  from  SMR  to  CMR  and  from  CMR  to  DMR, 
respectively,  are  taken  as  M2t  >  1  and  M2k  >  l. 

It  is  now  also  quite  certain  that  the  v on  Neumann 
paradox  can  be  explained  on  the  basis  of  the  viscous 
boundary-layer-displacement  thickness  at  the  reflec¬ 
tion  point  P  of  RR. 

Some  useful  discussions  on  the  effect  of  Mach- 
stem  curvature  and  its  influence  on  calculating 
precise  values  of  the  triple-point-trajectory  angle 
are  given. 
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where  i  and  j  are  the  upstream  and  downstream 
values,  respectively. 

In  the  case  of  regular  reflection,  or  the 
two-jhook  theory,  two  sets  of  equations  for  i  «  0, 
j  =  1  and  i  •  1,  j  «  2  are  solved  under  the  bound¬ 
ary  condition  6i  =  62-  This  means  that  the  direc¬ 
tion  of  the  flow  behind  the  reflected  shock  wave  is 
along  the  wedge  surface.  The  physical  validity  of 
the  boundary  condition  will  be  discussed  later  in 
Section  3.4.  In  the  case  of  Mach  reflection,  or 
the  threeshook  theory,  three  sets  of  equations 
i  ■  0,  j  *  1;  i  «  1,  j  •  2j  and  i  «  0,  j  »  3  (<ti 
is  called  $3  in  this  case)  are  solved  under  the 
boundary  conditions  63  *  0\ -62  and  P2  *  P3,  which 
mean  that  the  flows  are  parallel  and  the  pressures 
are  identical  on  both  sides  of  the  slipstream.  In 
the  case  of  MR,  it  is  also  assumed  that  the  Mach 
stem  is  perpendicular  to  the  wedge  surface,  which 
is  reasonable  but  not  precise  experimentally.  This 
assumption  will  be  discussed  in  Section  3.2.  The 
boundary  conditions  63  =  91-02  is  also  discussed  in 
Appendix  C.  The  computer  program  used  to  solve 
these  sets  of  equations  for  obtaining  the  transi¬ 
tion  boundaries  in  the  (Ms- 6^)  or  (Ms-6*)  planes 
is  shown  in  Ref.  30. 


Comments  are  also  presented  in  the  appendices 
on  the  nature  of  the  solutions  of  the  three-shock 
theory,  the  effect  of  slipstream  thickness  and  the 
transition  CMR  *  DMR. 


2.  ANALYSES 

2.1  Method  of  Calculations 

The  basic  equations  and  assumptions  used  in  the 
present  calculations  are  the  same  as  those  used  by 
Ben-Dor  (Ref.  3)  unless,  in  some  cases,  the  assump¬ 
tions  and  boundary  conditions  are  changed  to  evalu¬ 
ate  their  effects  on  the  solutions.  The  method 
emp,'  'yed  for  the  calculations  is  briefly  described 
below. 

The  phenomena  are  assumed  to  be  pseudo -stationary 
and  all  the  velocities  in  the  equations  below  are 
those  relative  to  a  reference  frame  attached  to  a 
point  which  moves  with  a  constant  velocity.  The 
reference  point  is  the  reflection  point  P,  in  RR, 
and  the  triple  point  T,  in  MR,  respectively  (see 
Fig.  2).  Each  region  which  is  divided  by  shock 
waves  and  a  slipstream  is  designated  by  0  to  3,  as 
shown  in  Fig.  2.  The  physical  quantities  on  both 
sides  of  each  shock  wave  in  the  vicinity  of  the 
reference  point  satisfy  the  following  equations: 

Conservation  of  tangential  velocity: 

oltan«l  «  p^tan^. -6.)  (1) 


2.2  Role  of  Relaxation  Lengths 

When  a  shock  wave  is  propagated  through  a  gas, 
the  translational  degrees  and  rotational  degrees 
of  freedom  of  the  gas  molecules  are  excited  to  the 
new  state  of  equilibrium  within  the  length  of  a 
few  mean-free-paths,  which  is  the  thickness  of  a 
shock  wave.  The  other  internal  degrees  of  freedom, 
however,  take  a  longer  time  to  reach  equilibrium. 

In  analyses  of  gasdynamic  phenomena  including 
shock  waves,  the  role  of  relaxation  lengths  is  very 
important.  That  is,  if  the  relaxation  length  of  an 
internal  degree  of  freedom  is  much  longer  than  a 
characteristic  length  of  the  phenomenon,  the  in¬ 
ternal  degree  of  freedom  can  be  treated  as  frozen 
at  the  initial  state.  If  an  internal  degree  of 
freedom  whose  relaxation  length  is  considerably 
shorter  than  a  characteristic  length  of  the  pheno¬ 
menon  it  can  be  assumed  to  be  in  equilibrium 
immediately  behind  the  shock  wave.  The  gas  is  in 
nonequilibrium  when  the  relaxation  length  and  the 
characteristic  length  lie  between  the  two  extreme 
cases  of  frozen  and  equilibrium  flow. 

in  the  problem  of  oblique -shock -wave  reflec¬ 
tions,  the  flow  Mach  number  at  the  second  triple 
point  is  used  as  a  criterion  for  the  transition 
from  CMR  to  DMR.  Therefore,  the  distance  between 
the  first  triple  point  T  and  the  second  triple 
point  T’  (Fig.  1)  can  be  considered  as  a  character¬ 
istic  length  since  the  state  at  the  second  triple 
point  depends  on  the  relaxation  process  which 
begins  at  the  incident  shock  wave.  Consider  the 
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angle  6  (Fig.  10),  which  is  the  angle  between  the 
incident  shock  wave  and  the  reflected  shock  wave. 

If  the  latter  is  curved  then  a  tangential  line  is 
drawn  at  the  triple  point  in  order  to  measure  S. 

In  this  case,  the  length  of  the  portion  of  the  re¬ 
flected  shock  wave  required  to  draw  a  tangential 
line  may  be  considered  as  another  characteristic 
length.  Define  a  practical  resolved  length  on  a 
photograph.  (Practical  means  that  it  is  different 
from  the  usual  resolution  determined  from  the 
ability  to  separate  two  close  points.)  The  prac¬ 
tical  resolved  length  depends  on  whether  the  photo¬ 
graph  is  an  interferogram,  schlieren  record  or  a 
shadowgram,  and  generally  it  is  harder  to  see  the 
exact  position  or  direction  of  a  shock  wave  on  an 
interferogram  than  on  a  schlieren  record  or  a 
shadowgram.  It  is  reasonable  to  choose  1  mm  as  a 
practical  resolved  length  on  an  interferogram  under 
the  conditions  of  the  experiments  in  the  present 
study. 

Strictly,  it  is  not  possible  to  determine  one 
single  characteristic  length  for  all  phenomena  of 
shock-wave  reflections.  The  measurement  of  an  angle 
has  its  characteristic  length,  which  may  be  the 
shortest  one,  and  the  criterion  for  transition 
between  CMR  and  DMR  has  another  characteristic 
length,  which  may  be  the  longest  one  in  the  present 
discussion.  The  two  characteristic  lengths  differ 
by  a  factor  of  about  10  in  a  typical  case.  If  the 
relaxation  length  has  an  intermediate  value  between 
these  two  characteristic  lengths,  the  phenomena 
cannot  be  analysed  by  the  simplified  method  used  in 
the  present  report,  which  assumes  the  internal 
degrees  of  freedom  as  being  either  in  equilibrium 
or  frozen.  It  is  reasonable  to  assume  that  the 
solution  of  this  case  lies  between  the  two  extreme 
cases. 

The  characteristic  lengths  have  an  uncertainty 
factor  of  about  10.  It  should  be  recalled  that  the 
relaxation  process  itself  does  not  have  a  precise 
length  for  it  takes  a  considerably  longer  distance 
than  the  relaxation  length  [defined  as  (1  -  1/e)  of 
the  final  value]  to  reach  the  final  equilibrium 
state.  This  does  not  mean  that  an  analysis  based 
on  relaxation  length  is  uncertain.  As  discussed 
later  in  this  section,  the  relaxation  lengths  vary 
easily  by  a  factor  of  103  'v  107  between  Ms  «  2  and 
Ms  =  10  under  the  same  initial  conditions,  which  is 
much  greater  than  the  uncertainty  of  the  character¬ 
istic  length  and  the  relaxation  process.  In  the 
discussions  which  follow,  a  characteristic  length 
of  1  mm  is  chosen.  This  is  essential  as  the  most 
quantitative  discussion  in  the  present  report  is 
based  on  the  values  of  various  angles.  Since  the 
relaxation  process  has  a  finite  length,  it  is  true 
that  the  angle  must  be  that  of  the  frozen-gas  case 
in  the  infinitesimal  vicinity  of  the  intersection 
point  of  two  shock  waves.  It  is  not  possible, 
however,  to  measure  the  angle  so  close  to  the 
intersection  point.  In  an  actual  measurement,  at 
least  1  mm  of  a  shock  wave  is  required  to  measure 
an  angle.  Consequently,  the  6  and  other  angles 
measured  are  those  between  a  tangential  line  in  a 
1  mm  region  when  the  shock  wave  is  curved.  Con¬ 
sequently,  whether  or  not  an  actual  transition  of  an 
internal  degree  of  freedom  from  the  frozen  to  the 
equilibrium  state  takes  place  in  some  finite 
region  over  the  transition  Mach  number,  is  deter¬ 
mined  from  the  condition,  whether  or  not  its  relaxa¬ 
tion  length  is  longer  than  this  characteristic 
length . 


In  the  present  report,  the  frozen  and  equili¬ 
brium-gas  assumptions  are  defined  as  follows.  The 
frozen-gas  assumption  means  that  only  the  transla¬ 
tional  and  rotational  degrees  of  freedom  are 
excited  to  their  equilibrium  values  and  that  the 
internal  degrees  of  freedom  are  all  frozen  at  their 
initial  states.  That  is,  y  *  1.667  for  Ar  and 
y  =  1.40  for  N2 ,  02  and  also  for  C02  (see  Appendix 
A).  As  discussed  later  in  this  section,  the 
internal  degrees  of  freedom  of  the  gases  considered, 
other  than  the  vibrational  degrees,  can  be  assumed 
as  frozen  at  their  initial  states  in  the  range  of 
the  experiments,  except  for  high  Mach  numbers  for 
C02 ,  since  the  relaxation  lengths  of  the  internal 
degrees  are  much  longer  than  the  defined  character¬ 
istic  length  of  the  phenomena.  Therefore,  in  this 
analysis  only  vibrational  excitation  is  the  addi¬ 
tional  internal  degree  of  freedom  required  for  the 
calculation  of  equilibrium-gas  properties,  as  com¬ 
pared  with  those  for  a  frozen  gas.  A  fictitious 
perfect-gas  case,  which  has  a  constant  y  of  1.29 
for  C02  is  also  computed,  even  though  it  has  no 
physical  significance  (see  Appendix  A) .  Neverthe¬ 
less,  the  agreement  of  the  C02  data  (experimental 
and  numerical)  in  the  (MS-6W)  or  (MS-9W)  planes  is 
excellent  by  accident  and  indicates  the  need  for 
better  criteria  for  SMR  CMR  and  CMR  *  DMR  tran¬ 
sition  lines. 

Figure  4  shows  the  vibrational  relaxation 
length  of  C02  (Ref.  17)  in  region  1  (the  region 
behind  the  incident  shock  wave)  at  the  initial 
conditions  p0  -  15  torr  and  Tq  =  300  K.  The 
pressure,  temperature  and  flow  velocity  in  region 
1  which  are  based  on  the  frozen-gas  assumption 
(>  »  1.40)  are  used  in  the  calculation  of  the 
relaxation  lengths.  The  relaxation  length  is  1  mm 
at  Ms  =  3.  Therefore,  C02  is  considered  to  change 
its  behaviour  from  a  frozen  gas  to  an  equilibrium 
gas  in  a  certain  region  at  Ms  =  3 .  Ando  (Ref.  14) 
discussed  the  vibrational  relaxation  of  C02  in  his 
paper.  He  misinterpreted,  however,  the  pressure 
behind  a  shock  wave  as  an  initial  pressure  when 
referring  to  the  experimental  results  of  Gaydon 
and  Hurle,  for  example.  Consequently,  he  obtained 
extremely  long  vibrational  relaxation  lengths  and 
concluded  that  the  vibrational  degrees  of  freedom 
of  C02  are  frozen  at  their  initial  states  under 
the  conditions  of  his  experiments. 

Figure  5  shows  an  interferogram  for  C02  at 
Ms  =  2.04,  po  ■  50  torr  and  To  =  297.3  K  taken  by 
Ando  (Ref.  14).  A  relaxation  process  behind  the 
incident  shock  wave  can  be  seen  clearly.  The 
vibrational  relaxation  length  at  Ms  =  2.04  is 
predicted  to  be  about  4  mm  according  to  Fig.  4. 

If  the  difference  in  the  initial  pressures  between 
the  case  shown  in  the  interferogram  and  in  Fig.  4 
(which  are  50  torr  and  15  torr,  respectively)  is 
taken  into  account,  the  relaxation  length  in  the 
case  of  the  interferogram  is  predicted  to  be  about 
1.3  mm.  The  relaxation  length  behind  the  incident 
shock  wave  in  Fig.  5  is  measured  as  1.7  mm,  which 
is  in  quite  good  agreement  with  the  prediction. 

The  dissociational  relaxation  length  of  C02  for 
the  same  initial  conditions  as  Fig.  4  is  about  70 
times  longer  than  its  vibrational  relaxation  length 
and  is  about  2  mm  at  Ms  =  10  (Ref.  17).  There 
might  be  some  effect  of  dissociational  excitation 
near  Ms  *  10  on  the  transition  boundaries  and  these 
would  have  to  be  included  for  given  initial  condi¬ 
tions  in  calculating  the  boundaries. 
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Figure  6  shows  the  vibrational  relaxation  length 
of  N2  (Refs.  18,  19).  The  initial  conditions  are 
the  same  as  Fig.  4  and  the  calculation  of  pj,  Tj 
and  u^  are  based  on  a  frozen-gas  assumption  Cy  =  1.4). 
Although  there  are  two  experimental  curves,  it  is 
probable  that  at  Mach  numbers  Ms  n,  io  the  vibra¬ 
tional  excitation  affects  the  flow  phenomena. 

Figure  7  shows  an  interferogram  of  N2  taken  by 
Ben-Dor  (Ref.  4)  at  Ms  =  8.06,  which  is  the  highest 
Mach  number  in  his  N2  experiments.  The  initial 
pressure  is  5.1  torr  and  temperature  is  298.2  K. 

A  relaxation  process  is  seen  behind  the  incident 
shock  wave.  (It  is  easier  to  see  if  a  tangent  line 
is  drawn  along  a  fringe.) 

The  dissociational  relaxation  length  in  the 
case  of  N2  is  much  longer  than  its  vibrational  re¬ 
laxation  length  like  ether  gases.  The  dissociation¬ 
al  relaxation  length  under  the  initial  conditions 
of  Fig.  6  is  about  200  mm,  even  at  Ms  »  15  (Ref.  20). 
Therefore,  the  dissociation  of  N2  can  be  neglected 
in  the  range  of  experiments  referred  to  in  the 
present  report . 

Figure  8  shows  the  vibrational  relaxation  length 
of  02  (Refs.  18,  19).  The  initial  conditions  and 
the  calculation  of  pj ,  Tj,  U[  are  the  same  as  those 
in  the  case  of  N2 .  The  vibrational  relaxation 
length  is  1  mm  at  about  Ms  =  6  n.  6.5  from  both 
references.  The  experiments  referred  to  in  the 
present  report  are  not  for  pure  02  but  air.  The 
vibrational  relaxation  length  of  02  in  air  may  have 
different  values  from  those  shown  in  Fig.  8,  as  the 
collisions  with  N2  molecules  are  more  frequent  than 
with  02 .  However,  it  is  reasonable  to  use  the  value 
for  pure  02  since  N2  is  considered  to  have  similar 
collision  properties  to  02 . 

Figure  9  shows  an  interferogram  of  a  double-Mach 
reflection  in  air  at  Ms  =  8.7,  pp  =  30.75  torr  and 
Tp  »  299  K  taken  by  Deschambault  (Ref.  21).  It 
shows  a  clear  relaxation  process  behind  the  incident 
shock  wave.  The  relaxation  length  is  not  measured 
since  it  is  difficult  to  determine  the  precise  point 
in  a  finite  interferogram  at  which  the  density  has 
approached  the  equilibrium  value  by  1/e  from  the 
frozen  value.  The  relaxation  process  seen  in  this 
interferogram  is  interpreted  as  that  of  N2  because 
the  relaxation  length  of  N2  at  the  conditions  of 
this  interferogram  is  predicted  to  be  about  2  mm, 
from  Fig.  6,  which  is  in  agreement  with  the  relaxa¬ 
tion  length  seen  in  the  interferogram.  No  inter¬ 
ferogram  was  available  showing  the  relaxation 
process  in  02 ,  which  occurs  at  a  lower  Mach  number 
than  that  in  N2 .  This  is  considered  due  to  the  fact 
that  the  concentration  of  02  in  air  is  only  20%,  so 
that  the  density  change  due  to  the  relaxation  of  02 
is  probably  too  small  to  be  seen  in  an  interferogram. 

The  dissociational  relaxation  length  of  02  at 
Ms  *  10,  under  conditions  of  Fig.  8,  is  about  30  mm 
(Ref.  20).  Consequently,  the  dissociation  of  02 
also  can  be  neglected. 

In  the  case  of  Ar,  since  Ar  is  a  monatomic  gas, 
there  is  neither  vibrational  excitation  nor  dissoc¬ 
iation.  Therefore,  Ar  can  be  treated  as  a  frozen 
gas  unless  electronic  excitation  and  ionization 
begin  to  take  place.  In  the  range  of  experiments 
referred  to  in  the  present  report,  electronic  ex¬ 
citation  and  ionization  can  be  neglected  and  Ar  can 
be  treated  as  a  frozen  gas. 


A  summary  of  the  discussions  in  this  section 
and  Section  3.1  are  given  in  Table  1. 


3.  RESULTS  AND  DISCUSSIONS 

3 . 1  Comparison  of  Experimental  and  Calculated 

Results  for  5 

As  expected,  the  calculated  results  based  on 
the  frozen-gas  model  differ  from  the  equilibrium- 
gas  model.  By  comparing  an  experimental  result 
with  the  frozen  or  equilibrium-gas  model  for  a 
given  flow  property,  we  can  determine  which  model 
is  valid.  It  is  important  to  choose  a  basic 
quantity  for  comparison  in  order  to  avoid  the 
effects  of  other  assumptions  used  in  the  calcula¬ 
tions.  From  a  practical  point,  the  chosen  property 
should  be  measured  easily  and  accurately.  Needless 
to  say,  the  conclusions  from  a  comparison  of  models 
must  be  consistent  with  the  discussion  on  relaxa¬ 
tion  lengths. 

The  angle  6  was  chosen  for  this  purpose.  As 
shown  in  Fig.  10,  it  is  the  angle  between  the 
incident  and  reflected  shock  waves.  The  comparison 
of  *  is  essentially  equivalent  to  the  comparison 
of  which  is  the  angle  between  the  incident 
shock  wave  and  the  triple-point-trajectory  path  at 
angle  v  (Fig.  10) .  The  comparisons  are  identical 
if  an  experimental  result  is  compared  with  a  cal¬ 
culated  result  which  has  the  same  value  of 
(since  6  =  90°  +  0^  -  The  use  of  f  or 

obtained  experimentally  is  not  identical,  if 
compared  to  a  calculated  result  which  has  the  same 
value  of  0„  as  that  of  the  experiment,  as  done  in 
the  previous  papers  (since  :  =  90°  ♦  fw  +  ,  -  ,■'). 

A  knowledge  of  «  is  then  required  and  at  present 
there  is  no  precise  way  of  calculating  it.  This 
will  be  discussed  further.  A  comparison  of  f  has 
some  advantages  over  when  comparing  results, 
for  f  is  a  monotonic  function  of  %  at  fixed  Ms . 
However,  w'  has  a  minimum  point  at  a  certain  6W 
for  some  part  of  the  Mach  number  range  (Ref.  23). 
Also,  since  6  will  be  used  as  a  new  criterion  for 
transition  between  SMR  and  CMR  (Section  3.3),  it 
is  necessary  to  check  the  agreement  between 
measured  values  of  6  and  the  calculated  results. 

The  important  difference  between  the  present 
and  previous  comparisons  (Refs.  14,  23)  is  that, 
in  the  present  comparison,  the  experimental  values 
of  6  are  compared  with  calculated  results  for  the 
same  0^,,  instead  of  6W.  It  is  because  the  two 
incident  wave  angles  of  a  Mach  reflection,  i>p  and 
4>j,  together  with  Ms  show  that  the  wave  system 
can  be  derived  from  9/  and  respectively 
(4>0  =  90°  -  0*.  i}>3  =  90°  -  \)  .  In  the  latter,  the 
Mach  stem  is  assumed  perpendicular  to  the  wedge 
surface  at  the  triple  point,  which  is  not  valid 
experimentally,  as  described  in  Section  3.2. 
However,  0^  and  x  (or  and  $3)  are  not  inde¬ 
pendent  of  each  other.  The  three-shock  theory 
allows  us  to  choose  only  one  parameter  from  0^, 

X  and  0W  besides  Ms.  Therefore,  8*  should  be 
chosen  as  the  parameter  since  it  has  a  valid 
relation  with  $0,  while  \  or  6W  have  only  an 
assumed  Mach-stem  perpendicularity  relation  with 
$3  or  both  $3  and  respectively.  Therefore, 
comparison  with  a  calculation  having  the  same 
0W  as  the  experiment  is  not  suitable  until  an 
accurate  prediction  of  \  is  found.  Consequently, 
if  we  compare  the  results  of  an  experiment  with  a 
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calculation  having  the  same  value  of  (<w  means, 
effectively,  that  we  are  comparing  ‘  from  two 
different  trajectory  directions.  Comparisons 
using  both  ways  are  shown  in  Figs.  11  and  12,  for 
the  case  of  CO, . 

It  is  important  to  note  that  the  experiments 
which  have  different  initial  pressures  have  differ¬ 
ent  relaxation  lengths  even  at  fixed  Ms  and  To.  In 
addition,  factors  such  as  the  magnification  factor 
of  the  photograph,  the  direction  and  spacing  fringes 
in  finite-fringe  interferograms  and  a  residual 
initial  fringe  in  infinite-fringe  interferograms 
should  be  considered  in  the  reduction  of  experi¬ 
mental  data. 

Figure  11  shows  the  experimental  results  for 
in  C02  (Ref.  14)  plotted  in  the  (f-Ms)  plane.  The 
number  beside  each  experimental  point  shows  the 
value  of  0*  measured  by  Ando.  The  values  of  4  a,, 
determined  from  the  given  fiw  and  . The  solid 
lines  and  the  broken  lines  are  the  calculated 
for  frozen  (.  =  1.40)  and  equilibrium  C02 ,  respect¬ 
ively.  Also  a  perfect- gas  case  (,  =  1.29)  is  shown 
by  dotted  lines.  The  experimental  results  show 
quite  good  agreement  with  the  calculated  results 
based  on  the  equilibrium-gas  assumption  for  the 
entire  Mach  number  and  0W  ranges  in  which  experi¬ 
mental  data  were  obtained.  At  low  Mach  numbers, 
the  difference  between  the  calculated  results  of 
each  gas  model  becomes  so  small  that  it  is  of  the 
same  order  as  the  experimental  error.  This  result 
is  consistent  with  the  conclusion  in  Section  2.2, 
that  C02  begins  to  change  its  behaviour  from  a 
frozen  gas  (,  =  1.40)  to  an  equilibrium  gas  at 
about  Ms  =  3  for  T(J  =  300  K. 

Gvozdeva  et  al  (Ref.  16)  compared  their  experi¬ 
mental  results  of  CO,  with  their  analysis  in  the 

plane  (according  to  their  notation,  . and 
■j°  are  equal  to  our  . 1  and  To.  respectively).  They 
concluded  that  the  imperfect  gas,  including  vibra¬ 
tional  excitation,  showed  good  agreement.  It  is 
consistent  with  the  present  result  although  their 
experiments  were  done  at  only  one  Mach  number  (they 
quote  a  shock-wave  velocity  of  1600  m/s). 

Figure  12  shows  a  plot  using  Mw  rather  than  rw 
for  the  equilibrium  and  per feet -gas  models.  The 
frozen-gas  case  (,  =  1.40)  was  omitted  as  it  gave 
results  in  poor  agreement  with  experiments.  The 
agreement  is  not  as  good  as  obtained  in  Fig.  11. 

This  illustrates  the  additional  errors  caused  by 
when  using  ‘w  rather  than  0^. 

Ando  (Ref.  14)  compared  his  experimental  data 
(replotted  on  Fig.  11)  with  several  models  ranging 
from  a  fictitious  perfect  gas  (t  =  1.29)  to  complete 
equilibrium  in  vibration  and  dissociation  for  .'. 
Since  he  compared  his  experimental  data  with  his 
calculated  results  for  the  same  flw,  rather  than  0*, 
the  agreement  was  not  very  good.  Instead  he  chose 
to  base  his  agreement  in  the  (Ms-nw)  or  (Ms-Qw) 
plane  transition  boundaries  as  indicators  of  agree¬ 
ment  with  analysis.  In  those  planes  the  agreement 
is  accidentally  good  and  he  concluded  that  the 
perfect-gas  model  (,  =  1.29)  was  the  best.  However, 
his  comparison  of  with  imperfect-gas  cases 
including  excitation  of  three  or  four  vibrational 
modes  showed  the  best  agreement.  This  contradiction 
was  not  resolved  in  his  work.  Consequently,  we 
must  conclude  that  the  SMR  *  CMR  and  CMR  *  DMR  tran¬ 
sition  lines  must  be  replaced  by  new  criteria  such 


that  better  agreement  is  obtained  in  the  (Ms-6w)  or 
(Ms-Uw)  planes  for  all  tested  gases. 

Figure  13  shows  a  similar  plot  for  N2  of 
Ben-Dor's  data  (Ref.  4)  in  the  (Ms-‘)  plane.  The 
number  beside  each  experimental  point  shows  the 
value  of  0*  measured  by  Ben-Dor.  In  this  case 
(and  for  Ar,  Fig.  15)  the  values  of  -  had  to  be 
measured  from  the  photographs  in  his  report,  since 
neither  the  values  of  nor  .'  were  listed.  The 
calculated  results  are  shown  for  a  frozen  gas 
(T  =  1.4)  (solid  lines)  and  a  vibrat ional ly  equi¬ 
librium  gas  (broken  lines)  for  each  fixed  ‘w. 

Here  the  agreement  is  not  as  clear-cut  as  in  C02  . 
The  experimental  joints  lie  closer  to  the  equili¬ 
brium  lines  for  -w  30°  and  closer  to  the  frozen 
N2  for  30°. 

Ando  (Ref.  22)  remeasured  the  values  of  for 
N2  from  Ben-Dor's  original  photographs  and  obtained 
values  different  from  those  shown  in  Ben-Dor's 
paper.  All  values  but  one  obtained  by  Ando  for  ^ 
have  the  same  values,  or  greater,  than  those  by 
Ben-Dor.  The  greatest  difference  was  1.5°.  Ando's 
measurements  are  shown  separately  in  Fig.  14.  The 
experimental  results  in  Fig.  14  exhibit  the  same 
general  features  as  in  Fig.  13.  It  is  worth  noting 
that  the  error  in  the  measurements  is  of  the  same 
order  as  the  difference  between  the  frozen  and 
equilibrium-gas  cases.  Consequently,  a  clear  con¬ 
clusion  cannot  be  drawn  from  this  comparison. 
However,  the  results  do  not  contradict  the  previous 
discussions  on  relaxation  lengths.  Therefore,  N2 
can  be  treated  as  a  frozen  gas  (y  =  1.40)  except 
at  high  Mach  numbers,  Ms  10.  However,  in  the 
case  of  C02  the  error  in  the  measurements  does  not 
play  as  important  a  role  since  the  differences 
between  the  frozen  and  equilibrium-gas  calculations 
are  larger  than  the  experimental  error  and  points 
to  the  validity  of  C02  being  in  equilibrium  for 
Ms  3. 

Figure  15  shows  the  (Ms--1)  plot  for  Ar  derived 
from  Ben-Dor's  data  (Ref.  4).  As  Ar  is  not  ex¬ 
cited  over  the  given  Ms  range  only  the  frozen-gas 
(t  *  1.667)  results  are  shown  as  solid  lines  for 
each  fixed  ’'w.  Good  agreement  is  obtained  with 
experimental  results,  although  the  experimental 
l's  are  slightly  larger  than  the  calculated  re¬ 
sults  at  higher  6W. 

Figure  16  shows  the  (Ms-r)  diagram  based  on 
the  experiments  of  Deschambault  (Ref.  21),  in  air. 
The  solid  and  broken  lines  are  the  calculated 
results  for  fixed  bK  for  air  as  a  frozen  ()  =  1.4) 
and  equilibrium  gas,  respectively.  The  experi¬ 
mental  results  do  not  show  any  clear-cut  picture. 
For  example,  at  rv  *  25°,  the  agreement  favours 
the  equilibrium  solution  up  to  Ms  ~  9.  Whereas, 
at  =  35°,  agreement  is  best  with  the  frozen 
line  up  to  Ms  9 .  The  equilibrium  lines  are 
also  favoured  for  4  <  Ms  <  6  at  =  35°  and  45°. 
Consequently,  the  experimental  results  for  air 
again  do  not  show  as  a  clear-cut  experimental 
decision  between  the  frozen  and  equilibrium-gas 
cases  in  N2  and  air.  Quantitatively  accurate 
measurements  are  required  since  the  differences 
between  the  two  models  are  not  as  large  as  in  C02 . 

It  should  be  noted  that  the  angle  4  was  cal¬ 
culated  on  the  basis  of  conditions  in  state  1. 
However,  state  2  is  also  involved.  Conditions  in 
state  2  depend  on  the  wedge  angle  snd  it  is 
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possible,  for  example,  for  state  1  to  be  frozen  and 
for  states  2  and  3  to  be  in  equilibrium.  Conse¬ 
quently,  the  present  computer  program  will  have  to 
be  modified  and  the  (Ms-5)  plots  will  have  to  be 
corrected  at  a  future  date.  Better  agreement  may 
result . 

3.2  Effect  of  Variation  of  Mach-Stern  Angle  at 

Triple  Point  from  being  Perpendicular  to 

Wedge  Surface 

The  available  analytical  results  for  the  triple¬ 
point-trajectory  angle  x  as  a  function  of  6W  and  Ms 
are  not  sufficiently  accurate.  For  example,  Ando 
(Ref.  14)  compared  his  predictions  of  x  with  the 
experimental  results  [see  Appendix  B,  Fig.  B1  (a)'v  (e) 
in  Ref.  14} .  The  experimental  data  did  not  agree 
with  his  predictions  based  on  various  models  of  C02 
excitation.  The  fictitious  perfect-gas  model 
(>  =  1.29)  accidentally  agreed  best  with  the  experi¬ 
mental  data,  and  is  inconsistent  with  the  concepts 
of  relaxation  lengths.  This  inconsistency  and 
fortuitous  agreement  can  be  partly  attributed  to 
the  assumption  that  the  Mach  stem  is  perpendicular 
to  the  wedge  surface  as  described  below. 

In  the  three-shock  theory,  the  effective  wedge 
angle  'J,  plays  an  important  role.  That  is,  once  a 
value  of  <}y  is  given  at  a  fixed  Ms,  the  directions 
of  the  reflected  shock  wave  and  the  Mach  stem  are 
determined  along  with  the  flow  variables  from  the 
solution  of  the  three-shock  theory.  In  this  se¬ 
quence,  a  value  of  the  actual  wedge  angle  0«  is 
obtained  if  the  Mach  stem  at  the  triple  point  is 
assumed  perpendicular  to  the  wedge  surface.  In 
other  words,  assumptions  other  than  the  perpendic¬ 
ular  condition  would  give  other  values  for  0W. 

This  means  that  is  a  function  of  4^  including 
the  assumption  concerning  the  orientation  of  the 
Mach  stem  to  the  wedge  surface.  Therefore,  the 
value  of  x  is  also  affected  by  the  orientation  of 
the  Mach  stem.  In  all  previous  calculations  of 
the  three-shock  theory  (Refs.  1,  2,  13,  23),  the 
Mach  stem  was  assumed  perpendicular  to  the  wedge 
surface  although  discrepancies  were  observed  ex¬ 
perimentally  owing  to  the  lack  of  a  better  one. 

This  assumption  is  satisfied  exactly  at  the  foot 
of  the  Mach  stem  to  meet  the  boundary  condition 
that  flows  on  both  sides  of  the  Mach  stem  must  be 
along  the  wedge  surface  (neglecting  the  effects  of 
the  viscous  boundary  layer).  It  would  also  be 
valid  at  the  triple  point  if  the  Mach  stem  were 
perfectly  straight. 

Experimentally,  however,  the  Mach  stem  gener¬ 
ally  is  curved  and  the  angle  f  between  the  Mach 
stem  at  the  triple  point  and  the  line  perpendicular 
to  the  wedge  (see  Fig.  10)  has  been  found  from 
present  experiments  to  have  a  concave  and  convex 
distribution  between  -3°  and  7.S0,  respectively. 
Figure  17  shows  an  example  of  a  Mach  stem  in  C02 
which  is  not  perpendicular  to  the  wedge  surface. 

The  angle  e  in  this  case  is  7.5°.  Calculations 
using  these  values  instead  of  r.«  0°  give  different 
results  for  The  results,  taking  into  account 
the  effects  of  e,  were  obtained  graphically  in 
order  to  save  computer  time.  Therefore,  the  re¬ 
sults  were  not  obtained  for  all  experimental  cases 
and  are  given  in  round  numbers  in  most  cases. 

Table  2  shows  the  measured  values  of  e  from  exist¬ 
ing  photographs  (Refs.  14,  4)  and  the  calculated 
values  of  r.  which  are  consistent  with  values  of 
both  0w  and  x  (and  consequently  with  flw).  The 
calculated  values  are  given  for  the  frozen,  perfect 


equilibrium-gas  cases.  The  experimental  results 
for  0W  less  than  20°  are  not  compared  with  calcula¬ 
tions  since  for  small  wedge  angles  an  error  of 
0.5°  in  x  corresponds  to  more  than  an  error  of 
3°  in  c.  Therefore,  the  accuracy  of  the  comparison 
is  poor. 

In  the  case  of  C02 ,  the  equilibrium-gas  case 
showed  better  agreement  with  experiments,  contrary 
to  the  comparison  based  on  the  assumption  of  c  =  0°. 
[Calculations  based  on  the  e  »  0°  assumption  agrees 
best  with  the  perfect  (y  *  1.29)  case.]  However, 
there  are  still  non-negl igible  discrepancies  bet¬ 
ween  the  calculated  and  experimental  values  of  c. 

In  the  case  of  N2 ,  the  equilibrium-gas  case  showed 
rather  good  agreement  with  experiments,  which  seems  » 

inconsistent  with  the  discussion  on  relaxation 
lengths.  In  the  case  of  Ar,  the  discrepancies  bet¬ 
ween  the  calculated  and  experimental  t  are  smaller 
than  those  of  N2  or  C02 .  This  is  due  to  the  fact  « 

that  c  is  close  to  0°  for  Ar.  The  discrepancies 
which  remain  even  if  the  effect  of  c  is  taken  into 
account  are  considered  to  arise  from  other  factors. 

For  example,  the  displacement  thickness  of  the 
slipstream  between  regions  2  and  3  violates  one  of 
the  boundary  conditions  of  the  three-shock  theory, 
and  is  discussed  in  Appendix  C.  Unfortunately,  it 
does  not  minimize  the  discrepancies. 

It  should  be  noted  that  the  value  of  c  has  a 
significant  effect  on  the  solutions,  especially  on 
,,  when  they  are  described  as  functions  of  Sw.  It 
means  that  discussions  arising  from  a  comparison 
of  experimental  results  with  solutions  of  the 
three-shock  theory  (which  includes  the  assumption 
that  ■  =  0°)  are  no  longer  accurate  if  they  are 
made  by  comparing  the  results  which  have  the  same 
value  of  w .  However,  the  assumption  that  e  =  0° 
does  not  affect  the  comparisons  if  the  experi¬ 
mental  results  are  compared  with  analysis  based 
on  the  same  value  of  0W.  As  noted  previously, 
c  affects  the  values  of  0W,  whereas  solutions 
based  on  9*  are  independent  of  t.  Therefore,  the 
transition  boundaries,  except  for  RR  *  MR,  are 
more  accurate  in  the  (Ms-”i)  plane  than  in  the 
(Ms-0w)  plane  as  long  as  the  available  predictions 
of  \  remain  imprecise.  However,  the  (Ms-?w)  plane 
suffers  from  the  limitations  that  the  transition 
boundary  lines  are  too  close  together,  especially 
at  higher  Ms,  and  it  is  more  difficult  to  differ¬ 
entiate  between  frozen  and  equilibrium-gas  bound¬ 
aries  as  well  as  to  determine  if  experiment  and 
analysis  of  the  various  regions  are  in  agreement. 

In  addition,  the  (Ms-0^)  plane  has  a  multi-valued 
portion  near  the  RR  *  MR  boundary  as  described  in 
Section  3.4.  A  precise  analytical  determination 
of  V  is  therefore  a  much  needed  future  require¬ 
ment  . 

3.3  A  New  Additional  Criterion  for  Transition 

from  SMR  to  CMR 

i 

In  the  previous  sections,  it  was  concluded  that 
C02  behaves  as  a  gas  in  equilibrium  except  at  low 
Mach  numbers.  As  shown  in  Fig.  18  there  are  sig¬ 
nificant  discrepancies  between  the  experimental  ’ 

data  of  the  three  types  of  MR  and  the  calculated 
transition  boundaries  based  on  an  equilibrium-gas 
assumption.  These  discrepancies  are  considered  as 
arising  from  inappropriate  transition  criteria. 

The  criteria  M2y  '  1  and  M2g  '  1  are  only  necessary 
conditions  for  the  transitions.  The  existence  of 
other  necessary  conditions  for  the  transition  from 
SMR  to  CMR  may  push  the  boundary  line  into  the  CMR 
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region  and  the  CMR  to  DMR  line  into  the  DMR  region. 
This  would  improve  the  agreement  significantly. 

Complex-Mach  reflection  has  been  characterized 
by  the  existence  of  a  band  of  compression  waves  at 
the  kink  behind  the  reflected  shock  wave  whether 
it  can  be  seen  in  an  interferogram  or  not  (Refs. 

10,  16).  The  condition  for  the  existence  of  a  band 
of  compression  waves  may  be  related  to  the  angle  f- 
between  the  incident  shock  wave  I  and  the  reflected 
shock  wave  R  at  the  triple  point  T  (Fig.  10),  as 
fol lows . 

Experimental  results  (Ref.  24)  have  shown  that 
the  ratio  between  the  distance  from  the  wedge 
corner  to  the  kink  along  the  horizontal  line  (  and 
the  distance  from  the  wedge  comer  to  the  incident 
shock  wave  along  the  horizontal  line  L  is  equal  to 
the  ratio  between  the  flow  velocity  behind  the 
incident  shock  wave  (in  the  laboratory  frame)  uj, 
and  the  incident  shock-wave  velocity  usi(,  except 
for  wedge  angles  near  the  RR  boundary.  That  is, 

■/L  =  uj,/us,  (see  Fig.  19a).  In  other  words,  the 
kink  is  moving  with  the  flow  velocity  in  region  1 
as  far  as  the  horizontal  component  is  concerned. 
Therefore,  it  is  seen  that  the  flow  velocity  in 
region  1  relative  to  the  kink  has  only  a  vertical 
component.  This  means  that  the  flow  moves  downward 
(or  parallel  to  the  incident  shock  wave)  as  seen 
from  a  frame  of  reference  attached  to  the  kink 
(Fig.  19b). 

If  we  call  the  incident  angles  and  the  densities 
behind  the  reflected  shock  wave  on  either  side  of 
the  kink  as  fig,  fjg,  .  . 2 .  as  illustrated  in  Fig. 

19c,  then  the  existence  of  a  band  of  compression 
waves  is  equivalent  to  the  relation  of  i  2  ani1 
is  also  equivalent  to  fjg-90°:  ■  !1ir-90°!  because 
the  flow  variables  upstream  on  either  side  of  the 
kink  are  the  same.  This  means  that  the  closer  to 
90°  to  the  shock  front  that  the  incident  streamlines 
are,  the  stronger  the  shock  wave,  and  consequently 
the  greater  the  density  change  across  the  wave. 

There  are  four  possible  geometrical  inclinations 
for  the  reflected  shock  wave  on  either  side  of  the 
kink  as  illustrated  in  Fig.  19d.  (It  must  be 
remembered  that  a  configuration  which  has  a  convex 
angle  to  the  upstream  flow  direction  cannot  maintain 
a  sharp  corner  like  a  kink.  Only  configurations  A 
and  B  of  the  four  satisfy  the  condition  Ujg-900 
I tlg-900  ,  and  in  both  A  and  B  the  angle  1  is 
greater  than  90°.  Therefore,  "  90°  is  considered 

a  necessary  condition  for  the  existence  of  a  band 
of  compression  waves,  which  is  required  for  the 
transition  to  CMR.  (It  is  only  a  necessary  condi¬ 
tion  because  configuration  C,  which  does  not  satisfy 
the  condition  fjg-go0,  •  ;$ig-90°!,  also  has  f  90°.) 

Figure  20  shows  a  plot  of  the  experimental 
values  of  *  with  shock  Mach  number  Ms  for  the 
various  types  of  reflections  for  N2,  Ar,  C02  and 
air.  All  experimental  CMR  but  one  in  C02  lie 
below  the  '  =  90°  line,  which  means  that  they 
satisfy  the  condition  A  90°.  This  provides  good 
evidence  that  A  •  90°  is  a  necessary  condition.  It 
must  be  pointed  out  that  results  in  Fig.  20  do  not 
include  any  calculation  or  assumptions.  They  show 
only  the  experimental  relation  between  the  SMR  *  CMR 
boundary  line  and  A.  It  means  that  this  is  a  direct 
comparison  between  the  new  necessary  condition  and 
the  experimental  data. 

As  mentioned,  the  new  condition  for  transition 
6  90°  is  only  a  necessary  condition.  Therefore, 


the  condition  M2T  a  1  remains  along  with  the  new 
condition.  The  transition  takes  place  when  both 
of  these  two  conditions  are  satisfied.  Conse¬ 
quently,  the  term  new  criterion  quoted  below 
includes  both  necessary  conditions,  and  the  criter¬ 
ion  M2T  1  is  the  former  criterion. 

Figure  21  shows  a  transition  boundary  plot  for 
equilibrium  C02  in  the  ( Ms -6w)  plane  with  emphasis 
on  the  SMR  *  CMR  transition.  The  solid  line  shows 
the  new  criterion  and  the  broken  line  shows  the 
former  criterion  M2T  *  1.  One  of  the  two  SMR 
points  which  was  calculated  in  the  CMR  region  of 
the  former  criterion  now  lies  in  the  SMR  region  of 
the  new  criterion. 

Figure  22  shows  a  similar  transition  boundary 
plot  in  the  (Ms-H^)  plane  of  frozen  N2  and  air 
emphasizing  the  SMR  *  CMR  transition.  The  solid 
line  and  the  broken  line  show  the  transition 
boundaries  between  SMR  and  CMR  based  on  the  new 
criterion  and  the  former  criterion,  respectively. 
Most  of  the  SMR  points  which  were  previously  in 
the  CMR  region  now  fall  into  the  SMR  region  of  the 
new  criterion.  The  experimental  distribution  of 
SMR  and  CMR  show  very  good  agreement  with  the  new 
criterion  line  except  for  two  points  in  air  at 
Ms  8.  This  tendency  is  consistent  with  the 
discussion  on  vibrational  relaxation  lengths  for 
02  and  N2 .  Figure  23  shows  the  calculated  transi¬ 
tion  boundary  lines  of  the  new  criterion  for 
frozen  (solid  line)  and  equilibrium  air  (broken 
line).  It  can  be  seen  clearly  in  this  figure 
that  the  experimental  boundary  of  air  between 
SMR  and  CMR  approaches  the  equilibrium-air  line 
at  high  Mach  numbers. 

Figure  24  shows  the  calculated  lines  for  M2T  =  1 
and  .*  =  90°  for  Ar.  In  this  case,  the  (  =  90°  line 
lies  below  the  M^x  =  1  line  for  Ms  ■  6  and  slightly 
above  for  Ms  6.  The  addition  of  the  condition 
4  90°  makes  little  difference  since  the  region 

above  both  lines  satisfies  both  necessary  condi¬ 
tions  and  therefore  the  line  which  lies  above  is 
the  new  transition  line,  which  is  close  to  the 
M21  -  1  line  in  this  case. 

It  can  be  concluded  that  the  new  criterion, 
composed  of  conditions  M?t  1  and  A  ''  90°,  gives 
a  transition  boundary  line  between  SMR  and  CMR 
which  agrei  11  with  the  distribution  of  experi¬ 
mental  points  for  all  gases  tested  using  calcula¬ 
tions  consistent  with  6  and  relaxation  lengths. 

There  are  also  discrepancies  in  the  transition 
boundaries  between  CMR  and  DMR.  In  the  case  of 
C02 ,  for  example,  as  shown  in  Fig.  18,  the  equili¬ 
brium-gas  assumption  gives  a  CMR  DMR  transition 
boundary  line  which  is  not  in  good  agreement  with 
the  experimental  result  if  the  boundary  line  is 
based  on  the  criterion  M2K  <  1 .  Results  for  other 
gases  also  have  discrepancies.  Perhaps  an  addi¬ 
tional  necessary  criterion  is  required  in  the  case 
of  the  transition  from  CMR  to  DMR.  However,  no 
additional  criterion  (M2K  1  is  one  of  the  neces¬ 
sary  conditions)  has  been  found  so  far.  Some 
comments  on  this  problem  are  given  in  Appendix  P. 

It  should  also  be  noted  that  the  CMR  *  DMR 
transition  line  is  based  on  the  assumption  that 
the  distance  (L-C)  between  the  two  triple  points 
T  and  T'  (see  Fig.  19a)  remains  finite.  As  noted 
above,  this  is  a  reasonable  assumption  for  the 
gases  investigated  as  long  as  0W  s  40°  and  Ms  2  2- 
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That  is,  beyond  this  region  it  is  found  experiment¬ 
ally  that  the  second  triple  point  T'  approaches  and 
merges  with  the  first  triple  point  T  (Refs.  13,  15). 
Consequently,  the  CMR  *  DMR  line  approaches  and 
becomes  coincident  with  the  SMR  *  CMR  line  at  the 
RR  *  MR  boundary  at  point  P,  as  illustrated  in  Fig. 
25.  No  analytical  solution  has  been  found  for  this 
behaviour.  It  can  be  seen  that  an  experimentally 
drawn  line  in  this  region  results  in  much  better 
agreement  with  the  remaining  data  ir.  the  CMR  and 
DMR  region. 

3.4  Persistence  of  Regular  Reflection  (RR)  into  the 

Region  of  Mach  Reflection  (MR)  —  von  Neumann's 

Paradox 

Experimental  results  show  that  RR  persists 
beyond  the  RR  *  MR  boundary  of  the  detachment 
criterion.  Figure  26  shows  the  persistence  of 
RR  in  the  case  of  C02 .  In  this  figure,  the  experi¬ 
mental  points  are  plotted  in  the  (Ms-9w)  plane 
instead  of  the  (Ms-0w)  plane  since  it  has  a  multi¬ 
valued  portion  near  the  RR  *  MR  boundary.  [At  the 
RR  boundary  of  the  MR  region,  \  has  a  finite  value. 
Therefore,  0*  -  ®w  *  X  has  a  larger  value  than  6W 
at  the  boundary.  On  the  other  hand,  6W  is  equal 
to  dw  in  the  RR  region.  Consequently,  the  RR  t  MR 
boundary  corresponds  to  two  separate  lines  (0W  and 
♦  x)  in  the  (Ms-Ow)  plane.  A  point  between 
these  two  lines  in  the  (Ms-B^)  plane  corresponds 
to  two  physical  points,  one  in  RR  and  the  other 
in  MR.)  Therefore,  a  discussion  on  the  persistence 
of  RR  is  best  made  in  the  (Ms-9w)  plane,  unlike 
the  boundaries  in  the  MR  region  which  are  best 
described  in  the  (Ms-0w)  plane  (see  Section  3.2). 
This  inconsistency  can  be  removed  in  the  future  if 
an  accurate  prediction  of  \  can  be  found,  then  the 
(MS-6W)  plane  could  be  used  throughout.  As  can  be 
seen  in  Fig.  26,  RR  persists  down  to  0W  =  42°  at 
Ms  =  7.2.  The  persistence  depends  on  Ms  and  the 
data  points  are  well  below  the  boundary  of  the 
detachment  criterion.  However,  the  experimental 
boundary  is  not  unique.  In  2  <  Ms  <  3,  DMR  points 
were  obtained  at  0W  *  48°  in  Ref.  26.  In  Ref.  22, 

RR  points  were  again  obtained  at  a  lower  wedge 
angle  of  9w  =  47°.  This  fact  suggests  that  the 
extent  of  the  persistence  depends  on  factors  other 
than  Ms  and  9W.  Figure  27  shows  the  persistence  of 
RR  in  the  case  of  air.  No  unique  experimental 
boundary  exists  in  this  case  also,  despite  the  fact 
that  the  persistence  is  very  apparent. 

Hornung  et  al  (Ref.  11)  suggest  that  the  persis¬ 
tence  of  RR  could  be  explained  in  terms  of  the 
viscous  boundary  layer  which  develops  behind  the 
reflection  point  P  (Fig.  28)  on  the  wedge  surface. 

In  the  viscous  boundary  layer  produced  after  the 
passage  of  the  reflection  point,  the  flow  has 
progressively  lower  velocities  (and  is  zero  at  the 
wall)  than  the  flow  in  the  free  stream  in  a  labora¬ 
tory  frame.  However,  the  flow  in  the  boundary 
layer  has  higher  velocities  than  the  flow  outside 
in  a  reference  frame  attached  to  the  reflection 
point.  This  corresponds  to  the  fact  that  the  wedge 
surface  is  moving  with  a  velocity  equal  and  opposite 
to  the  incident  shock  speed  from  the  reflection 
point  P.  As  a  result,  the  displacement  thickness 
of  the  boundary  layer  has  a  negative  value  (see 
Fig.  28).  Therefore,  the  deflection  angle  through 
the  reflected  shock  wave  is  less  than  required 
without  the  boundary  layer.  The  displacement  angle 
9(1,  or  the  angle  between  the  flow  affected  by  the 
boundary  layer  and  the  actual  wedge  surface, 
enables  RR  to  persist  to  lower  wedge  angles  because 


the  point  at  which  the  deflected  angle  through  the 
reflected  shock  wave  increases  to  the  maximum  de¬ 
flection  angle  and  02m  is  now  lower  than  that 
without  the  displacement  angle.  It  should  be 
pointed  out  that  the  displacement  angle  is  not 
equal  to  the  difference  in  the  transition  wedge 
angle  caused  by  the  boundary-layer  effect.  In 
fact,  the  displaced  wedge  has  a  larger  inclination 
despite  the  fact  that  the  displacement  effect 
lowers  the  transition  wedge  angle  and  also  the 
wedge  angle  in  front  of  the  incident  shock  wave 
is  still  equal  to  the  actual  wedge  angle  6W.  Since 
the  reflected  shock  wave  is  of  the  weak  family  of 
the  two  possible  solutions,  the  angle  u>’  is  smaller 
than  the  one  without  a  boundary  layer.  Hornung 
et  al  (Ref.  11)  mentioned  that  the  opposite  was 
observed  experimentally.  However,  the  present 
analysis  of  previous  experiments  gives  a  smaller 
u)'  than  obtained  from  a  calculation  without  the 
boundary-layer  displacement  effect,  as  expected. 

The  angle  u>'  is  a  very  good  indicator  of  the 
boundary-layer  displacement.  Figure  29  shows  a 
(Ms-u')  plot  for  C02  experiments  (Refs.  14,  22) 
compared  with  calculations  using  several  displace¬ 
ment  angles.  There  are  no  analytical  solutions  of 
iv '  above  the  lines  corresponding  to  the  RR  ?  HR 
transition  boundaries  (for  each  0d)  .  The  experi¬ 
mental  ui'  are  smaller  than  those  for  0d  =  0°.  The 
differences  from  the  calculated  line  9d  =  0°  are 
larger  when  the  wedge  angle  is  lower,  in  other 
words,  closer  to  the  MR  region.  This  agrees  quali¬ 
tatively  with  the  analytical  fact  that  the  effect 
of  the  boundary-layer-displacement  angle  on  v' 
becomes  smaller  at  higher  wedge  angles.  The  experi¬ 
mental  results  agree  with  a  displacement  angle  of 
9d  -  "1°  -2°.  A  quantitative  discussion  of  this 

value  is  given  subsequently.  Figure  30  shows 
similar  experimental  results  for  air  and  N2 .  For 
simplicity,  only  experimental  points  of  9W  =  50° 
are  plotted  and  compared  with  calculations.  Ik  th, 
erase  of  air ,  the  comparison  is  complicated  because 
air  and  Ns  changes  from  a  frozen  tc  an  equilibrium 
gas  at  an  intermediate  Mach  number.  (The  reason 
why  the  frozen  line  0d  *  0°  does  not  exist  at  high 
Mach  number  is  that  Bw  -  50°  with  0d  =  0°  (frozen) 
is  in  the  MR  region.)  At  lower  Mach  numbers,  the 
experimental  results  agree  approximately  with  the 
calculation  for  a  frozen  gas  with  0d  =  -1°.  At 
higher  Ms,  the  tendency  is  for  the  points  to 
approach  the  equilibrium  lines.  Experimental 
points  for  other  0W  are  compared  with  the  9d  =  -1° 
case  for  frozen  and  equilibrium  air  in  Fig.  31. 

Other  Bw  show  a  behaviour  similar  to  the  case  of 
9W  -  50°.  Figure  32  shows  the  results  for  Ar.  The 
two  available  experimental  points  agree  with  the 
calculation  of  0d  -  "2°. 

For  a  complete  analysis,  the  displacement  angles 
must  be  based  on  boundary-layer  theory  (Ref.  29). 

It  is  not  possible  to  obtain  a  unique  displacement 
angle  even  if  the  initial  pressure  is  given.  The 
displacement  thickness  itself  can  be  calculated  and 
is  shown  in  Table  3  for  the  case  of  air.  The  dis¬ 
placement  thickness  is  proportional  to  >x,  which 
means  that  the  slope  of  the  displacement  thickness, 
which  is  proportional  to  1/r'x,  becomes  infinite  in 
the  limit  as  x  -*  0.  In  other  words,  arbitrary 
angles  can  be  obtained  by  choosing  a  distance  from 
the  reflection  point  in  which  the  mean  slope  is 
calculated.  Here  we  face  a  problem  similar  to  the 
choice  of  a  characteristic  length  in  Section  2.2. 

An  appropriate  distance  should  be  chosen  as  the 
distance  for  calculating  a  mean  slope.  It  is 
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reasonable  to  choose  1  mm  as  the  distance,  similar 
to  the  case  of  the  relaxation  length.  The  mean 
slopes  for  1  mm  are  given  in  Table  3  for  the  RR z MR 
boundary  in  air.  The  angles  lie  between  -1°  and 
-2°  and  are  in  good  agreement  with  the  present 
comparison . 

The  RR  *  MR  transition  boundary  lines  with 
several  boundary-layer-displacement  angles  Od  are 
compared  with  the  experimental  persistence  of  RR 
in  Figs.  26  and  27  for  C02  and  air,  respectively. 

In  the  case  of  C02 ,  a  displacement  angle  of  -1° 
changes  the  transition  boundary  by  about  0.5°.  All 
experimental  RR  are  above  the  equilibrium-gas  case 
of  dj  =  -2°.  At  about  Ms  =  2,  which  is  in  the 
region  where  C02  could  be  a  frozen  gas,  RR  persist 
below  6d  *  -4°  of  the  frozen-gas  case  and  they 
seem  to  correspond  to  about  Bj  =  -5°.  The  transi¬ 
tion  from  RR  to  MR  is  dominated  by  the  flow  at  the 
reflection  point.  Although  it  is  meaningless  to 
assume  that  the  flow  properties  in  the  infinite¬ 
simal  vicinity  of  the  reflection  point  determines 
the  reflection  phenomena,  the  region  which  domin¬ 
ates  the  phenomena  is  considered  to  be  much  less 
than  1  mm,  which  is  the  characteristic  length 
chosen  in  the  discussion  on  relaxation  lengths. 
Consequently,  the  Mach  number  at  which  vibrational- 
excitation  effects  begin  to  take  place  in  the 
transition  is  higher  than  Ms  =  3.  For  the  same 
reason,  the  displacement  angle  which  dominates  the 
transition  must  have  a  different  value,  which  is 
larger  than  that  for  a  change  in  the  angle  u'1. 

Since  the  size  of  the  region  which  dominates  the 
reflection  phenomena  is  not  known,  no  further 
quantitative  comparison  is  possible.  The  variation 
of  pressures  must  be  taken  into  account  in  this 
discussion,  even  in  the  frozen-gas  case,  since  the 
thickness  of  the  boundary  layer  depends  on  the 
pressure  or  Reynolds  number.  Figure  27  shows  the 
same  comparison  for  air.  Since  air  behaves  as  a 
frozen  gas  in  the  region  shown  in  Fig.  27,  only 
the  frozen-gas  cases  are  shown.  Regular  reflection 
(RR)  persists  down  to  the  transition  line  of  84 =  -2° 
at  lower  Mach  numbers  and  Bj  =  -4°  at  higher  Mach 
numbers.  A  further  quantitative  analysis  is  not 
possible  for  the  reasons  already  given  for  C02  . 

Detailed  calculations  of  the  flow  quantities 
in  the  various  states  generated  by  pseudo-stationary 
oblique-shock-wave  reflections  in  tabular  form,  for 
the  convenience  of  the  researcher,  can  be  found  in 
Ref.  30. 


4.  CONCLUSIONS 

The  recent  results  on  pseudo-stationary  oblique- 
shock-wave  reflections  were  assessed  based  on 
criteria  associated  with  molecular  relaxation 
lengths,  the  basic  angle  6  between  the  incident 
shock  wave  and  the  reflected  shock  wave  and 
transition  boundaries  in  the  (Ms-6^)  plane. 

The  following  conclusions  apply  in  the  range 
1  <  Ms  <  10;  at  Mach  numbers  lower  than  3,  C02 
behaves  as  a  frozen  gas  (y  *  1 .40)  and  as  an 
equilibrium  gas  with  vibrational  excitation  for 
Ms  >  3;  nitrogen  and  air  are  frozen  (y  »  1.40) 
below  M*  *  6  n,  9,  depending  on  the  initial  pressure; 
argon  is  frozen  in  the  present  Mach  number  range. 
Dissociation  and  ionization  can  be  neglected 
(except  for  C02  at  lower  pressures)  up  to  Mg  *  10. 


The  curvature  of  the  Mach  stem  affects  the 
triple-point-trajectory  angle  x  significantly. 

This  discrepancy  accounts  for  most  of  the  inaccu¬ 
racies  of  predicting  x-  There  are,  however,  still 
non-negl igible  discrepancies  between  the  experi¬ 
mentally-measured  x  and  the  calculated  x  even  if 
the  effect  of  Mach-stem  curvature  is  taken  into 
account . 

The  condition  that  the  angle  between  the 
incident  shock  wave  and  the  reflected  shock  wave 
is  greater  than  90°  is  an  additional  necessary 
condition  for  the  transition  from  SMR  to  CMR.  The 
addition  of  the  new  criterion  improves  the  agree¬ 
ment  between  analysis  and  experiments. 

The  persistence  of  RR  is  explained  in  terms  of 
the  induced  viscous  boundary  layer  on  the  wedge 
surface.  The  measured  reflected-wave  angles  w' 
agree  with  the  calculated  values  which  take  into 
account  the  boundary-layer-displacement  effect. 

The  value  of  the  boundary-layer-displacement  angle 
9d  necessary  to  explain  the  persistence  of  RR 
is  reasonable  compared  with  the  calculation  of 
u )' .  Rigorous  quantitative  predictions  are 
difficult  owing  to  the  uncertainty  in  calculating 
the  boundary- layer-displacement  angle. 

Large  discrepancies  at  the  CMR  *  DMR  boundary 
between  calculated  and  experimental  results  remain. 
In  addition,  the  experimental  results  show  that  the 
CMR  Z  DMR  1 ine  and  the  SMR  Z  CMR  1 ine  merge  at 
the  RR  line  due  to  the  merging  of  the  two  triple 
points  in  DMR.  An  analytical  prediction  for  this 
effect  is  not  available.  Analytical  predictions 
for  Mach-stem  curvature  and  its  relation  to  the 
triple-point-trajectory  angle  x  at  low  shock  Mach 
numbers  are  also  required. 
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Table  1 


Range  of  experimental  pressure  and  shock  Mach  number 
where  gases  can  be  treated  as  frozen  or  equilibrium 


for  oblique-shock -wave  reflections 


*at  pQ  *  IS  torr 
**at  pQ  =  SO  torr 


Table  2 


Comparison  of  calculated  and  experimental  values  of  e 


(experimental 

values 

are  measured 

from  Ref. 

Gas 

0 

w 

(deg) 

M 

s 

e 

(exp) 

(deg) 

Frozen 

(deg) 

e 

(cal) 

Perfect 

(deg) 

Equilibrium 

(deg) 

C02 

20.3 

4.72 

7.5 

-0.5 

1.5 

3.5 

20.3 

6.0 

6.5 

-1.0 

1.0 

3.5 

20.3 

10.18 

4.5 

-0.5 

1.5 

4.5 

27.1 

1.78 

-1.0 

-5.5 

-3.0 

-2.0 

27.1 

2. SO 

4.0 

-2.0 

0.0 

1.0 

27.1 

3.47 

5.5 

0.0 

1.5 

3.0 

27.1 

4.05 

4.0 

-O.S 

1.0 

2.5 

27.1 

8.16 

6.5 

-1.0 

1.0 

3.5 

30.1 

2.04 

2.0 

-3.0 

-1.0 

0.0 

30.1 

3.44 

4.5 

0.0 

1.5 

2.5 

30.1 

6.32 

6.0 

-0.5 

1.0 

3.0 

30.1 

8.30 

7.0 

0.0 

1.5 

3.5 

30.1 

9.77 

6.5 

-2.0 

0.0 

2.5 

N2 

20 

1.93 

0.0 

-5.0 

— 

-4.5 

20 

3.74 

2.0 

1.0 

- 

2.0 

20 

4.81 

S.O 

2.5 

- 

4.0 

20 

6.87 

4.0 

2.0 

- 

4.0 

30 

1.97 

0.0 

-1.0 

- 

1.0 

30 

3.68 

2.0 

1.5 

- 

2.0 

30 

S.93 

3.0 

1.5 

— 

3.0 

30 

7.97 

4.0 

0.5 

- 

2.0 

Ar 

20 

2.00 

-2.0 

-2.5 

_ 

20 

4.40 

O.S 

0.0 

— 

— 

20 

5.20 

0.0 

-0.5 

- 

— 

20 

7.76 

0.5 

0.0 

— 

— 

30 

2.03 

-2.S 

-2.0 

— 

— 

30 

4. SI 

0.5 

0.0 

— 

- 

30 

6.36 

0.0 

O.S 

- 

— 

30 

8.01 

1.0 

0.0 

— 

- 

40 

2. OS 

0.0 

1.0 

- 

— 

40 

S.28 

0.5 

-1.0 

- 

— 

40 

6.o2 

1.0 

-1.0 

— 

— 

Table  3 

Boundary-layer-displacement  thickness  6*  and  boundary-layer-displacement  angle  8 

on  wedge  surface 

(frozen  air  at  RR  t  MR  boundary,  pQ  =  15  torr,  Tq  =  300  K,  based  on  Ref.  29) 


M 

s 

f>*/v7 

cm'1 

(mean  value 
over  1  mm) 

1.2 

-5 . 83xlO~3 

-1. 

06 

1.4 

-6.74xl0"3 

-1 

22 

1.6 

-6.82xl0"3 

-1 

24 

1.8 

-6.94xl0"3 

-1. 

26 

2.0 

-6 ,90xl0"3 

-1. 

25 

2 .  S 

-6 . 79xl0"3 

-1 

23 

5.0 

-6 . 8Sxl0"3 

-1 

24 

4.0 

-7 . lOxlO3 

-o. 

29 

5.0 

-7 . SSxlO3 

-1 

37 

7.0 

-8.41xl03 

-1.52 

10.0 

-9.84xl0'3 

-1 

78 

FIG.  2  REGULAR  AND  MACH  REFLECTIONS  IN  STATIONARY  COORDINATES  WITH 
RESPECT  TO  POINTS  P  AND  T,  RESPECTIVELY. 

0,  1,  2,  3  -  REGIONS;  I  -  INCIDENT  SHOCK  WAVE;  R  -  REFLECTED 
SHOCK  WAVE;  M  -  MACH  STEM;  P  -  REFLECTION  POINT;  T  -  TRIPLE 
POINT;  8  -  DEFLECTION  ANGLE;  <f>  -  WAVE  ANGLE;  rmnnn  WEDGE; 

-*  -  FLOW  DIRECTION;  -  SLIPSTREAM;  -  -  DIRECTION  OF 

TRIPLE -POINT  PATH. 


e 


0 

} 

i 

FIG.  3  SHOCK-POLAR  ILLUSTRATING  RR  t  MR  TRANSITION  BASED  ON  THREE 
DIFFERENT  CRITERIA. 


I  -  INCIDENT  SHOCK  WAVE;  R,  R',  R"  -  REFLECTED  SHOCK  WAVE  AT 
RR  Z  MR  TRANSITION  (R  -  DETACHMENT  CRITERION  AT  A  -*■  B;  R'  - 
MECHANICAL-EQUILIBRIUM  CRITERION  AT  C;  R"  -  SONIC  CRITERION 
AT  D  -*•  E) ;  A,  B,  C,  D,  E  -  STATES  BEHIND  REFLECTED  SHOCK  WAVE 
AT  RR  t  MR  TRANSITION  (A  -  RR  STATE  OF  DETACHMENT  CRITERION; 

B  -  MR  STATE  OF  DETACHMENT  CRITERION;  C  -  RR  AND  MR  STATE  OF 
MECHANICAL-EQUILIBRIUM  CRITERION;  D  -  RR  STATE  OF  SONIC 
CRITERION;  E  -  MR  STATE  OF  SONIC  CRITERION). 


VIBRATIONAL  RELAXATION  LENGTH  (mm) 


FIG. 


SHOCK  MACH  NUMBER  Ms 


4  VIBRATIONAL  RELAXATION  LENGTHS  BEHIND  SHOCK  WAVES  IN  C02 . 

INITIAL  CONDITIONS:  p0  =  15  TORR,  T0  =  300  K  (REF.  17). 

FOR  HIGHER  OR  LOWER  INITIAL  PRESSURES  AT  300  K  THE  RELAXATION 
LENGTHS  ARE  PROPORTIONATELY  DECREASED  OR  INCREASED,  RESPECTIVELY. 
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FIG.  5  SINGI.H-MACII  REFLECTION  IN  C02  SHOWING  VIBRATIONAL  RELAXATION 

PROCESS  BEHIND  INCIDENT  SHOCK  WAVE  (FINITE-FRINGE  INTERFEROGRAM) . 
0W=  50°.  Ms  =  2.04,  p„  =  50  TORR,  T0  =  297.5  K  (REF.  14). 


ON  IN  N,  SHOWING  VIBRATIONAL  RELAXATION 
ENT  SHOCK  WAVE  (ORIENTATION  OF  FRINGES 
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VARIATION  of  ANr/F  *  .  S 

-°_=  3p0  K>  Po  =  IS  T0RRSH°CK'MACH  DUMBER  AT  FlXFn  a' 

*  F¥Dcf^UlLIBRIUM  CO,  ’  —  FROZEN  CO,  Cv  XE?  0*  F0R  CO;. 

experimental  aiW  ^ 


SHOCK  MACH  NUMBER  Ms 


VARIATION  OF  ANGLE  6  VS  SHOCK-MACH  NUMBER  AT  FIXED  0*  FOR  N2 . 

To  =  300  K,  po  =  15  TORR,  -  FROZEN  N2  (Y  =  1.4), 

-  EQUILIBRIUM  N2 ,  ♦  EXPERIMENTAL  DATA  FOR  6'  (0^  MEASURED 

BY  ANDO,  REF.  22) . 


*  / 


no 

100 -T47/  // 

/  /t36i  +32  6 

/  /  /X 

90 -J  +37^/  +32  f 
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EXAMPLE;  WHERE:  ENTIRE  MAGI  STEM  IS  NOT  PERPENDICULAR  TO  WEDGE 
SURFACE.  0W  =  20°,  Ms  =  4.72,  pD  =  14.5  TORR,  T0  =.297.7  K, 
CO2  (REF.  14) 


SHOCK  MACH  NUMOt 

0F  ™WS'T10N 
“&SISW  RESm-W  FOR  »sJw>  ,  .  CMR. 

^XPURlMF.NTS:  0 


bound arifs 
D  -  DMR  CM 


FIG.  19(a)  SCHEMATIC  DIAGRAM  ILLUSTRATING  ASSUMPTION  OF  POSITION  OF  KINK 
OR  SECOND  TRIPLE  POINT  AND  FIRST  TRIPLE  POINT. 


flow  vector  in  region  1 


FIG.  19(b)  SCHEMATIC  DIAGRAM  ILLUSTRATING  FLOW  DIRECTION  IN  REGION  1 
RELATIVE  TO  KINK. 


FIG.  19(d)  FOUR  POSSIBLE  GEOMETRIC  CONFIGURATIONS  OF  REFLECTED  SHOCK  WAVE 
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FIG.  22  SMR  CMR  TRANSITION  BOUNDARIES  OF  "NEW"  AND  "FORMER"  CRITERIA 
FOR  FROZEN  N2  AND  AIR  (y  =  1.40). 

-  NEW  CRITERION,  —.  —  —FORMER  CRITERION.  EXPERIMENTS: 

N2:  A  -  SMR,  V  -  CMR  (REF.  4);  AIR:  A  -  SMR,  ?  -  CMR  (REF.  21). 


I 

FIG.  23  COMPARISON  BETWEEN  FROZEN  AND  EQUILIBRIUM  AIR  OF  SMR  £  CMR 
I  TRANSITION  BOUNDARIES  OF  NEW  CRITERION. 

-  FROZEN  AIR,  —  EQUILIBRIUM  AIR.  EXPERIMENTS: 

A  -  SMR,  ▼  -  CMR  (REF.  21). 


SHOCK  MACH  NUMBER  Ms 


M2T  =  1  line  AND  6  =  90°  LINE  OF  FROZEN  ARGON 

— —  m2T  =  - 6  =  90°.  EXPERIMENTS:  A 

V  -  CMR  (REF.  4). 


SHOCK  MACH  NUMBER  Ms 


FIG.  25  EXPERIMENTAL  CONVERGENCE  OF  CMR  t  DMR  BOUNDARY  WITH  SMR  *  CMR 
BOUNDARY  AT  RR  *  MR  BOUNDARY  FOR  C02 . 

CALCULATED  BOUNDARY:  —  ■  —  FROZEN  C02  ,  —  ——EQUILIBRIUM  C02  . 
EXPERIMENTAL  POINTS:  o  -  RR,  A  -  SMR,  V  -  CMR,  □  -  DMR  (REF.  26) 

•  -  RR,  ▼  -  CMR,  ■-  DMR  (REF.  14) 

r  -  RR,  *  -  SMR,  tr~  DMR  (REF.  22) 

?  -  CMR,  a  -  DMR  (REF.  25) 

— - - EXPERIMENTAL  BOUNDARY  BETWEEN  CMR  AND  DMR. 
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FIG.  31  COMPARISON  OF  EXPERIMENTAL  REFLECTED  ANGLE  W  WITH  CALCULATIONS  OF  BOUNDARY- LAYER- 
DISPLACEMENT  ANGLE  6d  =  -1°  FOR  AIR  AND  N2 . 

-  FROZEN  AIR,  -  EQUILIBRIUM  AIR,  T0  =  300  K,  p0  =  15  TORR. 

EXPERIMENTAL  DATA:  AIR:  +  -  0W  =  45°,  ■  -  9W  =  47° ,  ▼  -  0W  =  49° , 

a  -  ew  =  50°,  •  -  ew  =  60°  (REF.  21), 


FIG.  32  COMPARISON  OF  EXPERIMENTAL  REFLECTED  WAVE  ANGLE  w'  WITH  CALCULATIONS 
WITH  SEVERAL  BOUNDARY-LAYER-DISPLACEMENT  ANGLES  0^  AT  0W  •=  C0°  FOR  Ar. 
o  -  EXPERIMENTAL  DATA  (REF.  4). 
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APPENDIX  A 


COMMENT  ON  •  rSr  FF  C02 


The  term  perfect  CO2  was  used  previously  in  Refs. 
13  and  14  to  indicate  a  gas  with  a  constant  value  of 
V  -  1.29  independent  of  temperature.  However,  from 
the  following  discussion,  CO2  with  a  constant 
1  ■  1.29  has  no  physical  validity.  At  room  tempera¬ 
ture,  the  vibrational  modes  of  CO2  are  partly  ex¬ 
cited.  The  doubly-degenerate  bending  modes,  which 
have  a  characteristic  temperature  of  960.2  K ,  are 
excited  to  14%  of  their  full  excitation  and  have  a 
contribution  of  0.9R  to  Cv-  The  other  two  stretch¬ 
ing  modes  contribute  O.tr  to  Cv.  Consequently, 
along  with  contributions  from  the  translational 
modes  and  rotational  modes  Cv  *  3.47R  at  room  tem¬ 
perature,  or  >  =  Cp/Cv.  =  1.29  at  room  temperature. 
However,  behind  a  shock  wave  the  translational  and 


rotational  temperatures  immediately  jump  to  the 
Kankine-Hugoniot  temperature  If.  Since  the  vib¬ 
rational  energy  is  too  slow  to  participate,  it 
remains  fro:en  on  either  side  of  the  shock  front 
and  the  temperature  T]f  can  be  computed  using 
1  =  1.4  |Ref.  31).  After  a  period  (relaxation 
time)  the  vibrational  modes  are  activated  and 
share  in  the  energy  through  equipartition.  The 
final  equilibrium  temperature  Tle  •  Tip.  :-.rf  •• 
CO2  with  •  *  1.29  has  no  physical  basis.  However, 
it  is  sometimes  used  in  engineering  applications 
to  give  an  approximate  answer  for  shock  waves  in 
CO2  with  vibrational  excitation.  The  error 
becomes  increasingly  worse  with  rising  shock 
Mach  number  (Ref.  31). 


APPENDIX  B 


THREE-SHOCK  THEORY 


The  three-shock  theory  provides  a  simple  and 
effective  analys.s  of  obi ique -shock -wave  reflec¬ 
tions.  Even  so,  it  includes  some  aspects  which 
are  physically  unreasonable  or  inconsistent  with 
experimental  results.  They  are  described  below. 

B.l  Behaviour  of  .  as  Ms  -  1 

The  values  of  v  at  fixed  9W  are  approximately 
constant  at  high  shock  Mach  numbers  Ms,  both 
experimentally  and  analytically.  As  Ms  decreases, 
the  values  of  ,  increase  except  for  large  ‘‘*w. 
Experimentally,  >  has  maximum  values  near  Ms  =  2 
and  decreases  for  other  Ms.  The  calculated  results 
of  <  from  the  three-shock  theory,  however,  do  not 
have  maxima  for  different  9W  but  continue  to  in¬ 
crease  as  Ms  ->  1,  as  shown  in  Fig.  A-l. 

Ben-Dor  (Ref.  3)  proposed  an  alternative  method 
to  obtain  x  at  small  wedge  Angles  (9*  <  4°)  which 
does  give  maxima  and  may  agree  with  experiment 
despite  his  failure  to  obtain  solutions  of  the 
three-shock  theory  for  small  0W  and  low  Ms .  His 
method  is  to  assume  that  0^  at  small  dw  coincides 
with  the  value^t  the  no  reflection  (NR)  boundary 
ywNR  an<*  X  =  9«NR  '  9w  at  small  wedge  angles.  Des¬ 
pite  that  his  method  gives  a  reasonable  approxima¬ 
tion  to  experiment  results  and  is  useful  from  an 
engineering  point  of  view,  he  gave  no  physical 
explanation  that  0^  at  small  0W  coincides  with  the 
NR  boundary.  His  method  does  not  give  the  neces¬ 
sary  dependency  of  °w  <=  f: wNR >  for  different  0W. 
Since  his  suggestion  is  outside  the  three-shock 
theory  realm,  being  only  an  assumption,  it  does 
not  provide  any  clue  for  improving  the  three-shock 
theory  or  an  alter:  !ve  theory  for  the  calculation 
of  /.  The  tendency  .0  increase  >  as  Ms  •»  I  is  not 
considered  realistic,  not  only  because  it  does  not 
agree  with  experimental  results,  but  it  also  yields 
a  no-solution  region  in  the  (MS-0W)  plane,  as  des¬ 
cribed  in  Appendix  B.2. 


The  reason  why  the  three- shock  theory  fails  to 
give  reasonable  solutions  at  low  Mach  numbers  ma\ 
be  as  follows.  The  reflected  shock  wave  is  com¬ 
posed  of  two  parts,  one  is  a  reflected  shock  wave 
emanating  from  the  triple  point  and  the  other  is 
a  bow  shock  wave  caused  by  the  interaction  of  the 
incident  flow  with  the  flow  deflected  by  the  wedge. 
The  two  shock  waves  intersect  at  the  kink  or  the 
second  triple  point  in  CMR  or  I'MR,  respectively. 

In  SMR,  the  two  shock  waves  blend  smoothly  so  that 
an  intersection  point  cannot  be  seen.  The  three- 
shock  theory  gives  the  position  of  only  one  of  the 
two  shocw  waves,  which  is  the  one  emanating  from 
the  triple  point.  Ihc  actual  shape  of  the  entire 
reflected  shock  wave  is  considered  to  be  determined 
by  the  two  shock  waves  mentioned  above  and  their 
interaction.  At  low  Mach  numbers,  both  the  reflect 
ed  shock  wave  from  the  triple  point  and  the  bow 
shock  wave  are  weak.  The  former,  however,  is  con¬ 
sidered  to  become  weaker  more  rapidly  than  the 
latter.  Therefore,  the  position  of  the  whole  re¬ 
flected  shock  wave  is  controlled  hy  the  bow  shock 
wave  rather  than  the  reflected  shock  wave  from 
the  triple  point.  In  the  limit  as  the  shock  Mach 
number  of  the  bow  shock  wave  approaches  1  (Mach 
wave),  or  w  *  0°,  the  bow  shock  will  intersect 
the  incident  shock  wave  at  a  point  =  "wNR>  w*1>c'h 
is  at  a  lower  (han  the  limit  of  solution  of  the 
three-shock  theory. 

The  reason  is  illustrated  in  Fig.  A -2.  Point 
0  in  Fig.  A-2  was  at  the  wedge  comer  when  the 
incident  shock  wave  collided  with  the  wedge  corner 
and  it  moves  with  the  flow  in  region  1_.  Therefore, 
the  flow  in  region  1  has  a  direction  TO  relative 
to  T  and  the  bow  shock  wave  in  the  vicinity  of  T  is 
normal  to  the  flow.  Since  the  bow  wave  is  a  Mach 
wave,  the  flow  velocity  in  region  1  relative  to  the 
bow  wave,  or  T,  is  sonic.  If  T  is  at  0^NR,  t*'e 
flow  Mach  number  in  region  1  is  also  1.0  according 
to  the  criterion  of  the  NR  boundary.  It  means  that 


T,  which  is  the  intersection  point  of  the  bow  shock 
wave  and  incident  shock  wave  in  the  limit  of  bow- 
shock  Mach  number  -»  1 ,  is  at  Therefore,  it 

is  reasonable  to  assume  that,  at  low  incident  Mach 
numbers  and  small  wedge  angles,  the  intersection  . 
point  of  the  extension  of  the  bow  shock  wave  and 
the  incident  shock  wave  is  nearer  a  point  for  which 
it  has  Consequently,  if  the  shape  of  the 

reflected  shock  wave  is  dominated  by  the  bow  shock 
wave,  the  bow  shock  wave  will  pull  down  the  posi¬ 
tion  of  the  reflected  shock  wave  emanating  from 
the  triple  point  and  eventually  the  position  of  the 
triple  point  itself  is  also  pulled  down  from  the 
position  given  by  the  solution  of  the  three-shock 
theory.  This  results  in  a  smaller  <  than  given  by 
the  solution  of  the  three-shock  theory. 

It  can  be  concluded  that  the  three-shock  theory 
is  inadequate  when  the  shock  wave  emanating  from 
the  triple  point  does  not  have  a  dominant  effect 
on  the  whole  shape  of  the  reflected  shock  wave.  In 
the  region  near  the  RR  t-  MR  boundary,  experimental 
results  show  that  the  position  of  the  kink  i/l  is 
no  longer  equal  to  the  ratio  between  the  flow  velo¬ 
city  behind  the  incident  shock  wave  and  the  velocity 
of  the  incident  shock  wave  ul;/us,.  This  is  also 
considered  to  be  caused  by  the  bow  shock  wave. 
Therefore,  calculations  which  take  into  account  the 
effects  of  the  bow  shock  wave  are  necessary  as 
Ms  -  1. 

B.2  No-Solution  Region  in  the  (M_.--^)  Plane 

The  regions  in  the  (Ms-0,i)  plane  have  an  addi¬ 
tional  region,  compared  with  those  in  the  (Ms-“w) 
plane,  namely  a  no-solution  (NR)  region.  The  flow 
Mach  number  behind  the  incident  shock  wave  relative 
to  the  triple  point  must  be  Mj  1  for  the  existence 
of  a  reflected  shock  wave  at  the  triple  point.  At 
a  fixed  Ms,  Mj  decreases  with  decreasing  0W  and  at 
a  certain  Mi  =  1 .  No  reflected  shock  wave  can 
exist  beyond  this  point.  The  line  Mi  =  1  is  the 
boundary  of  the  NR  region. 

it  was  thought  that  there  are  no  solutions  of 
the  three-shock  theory  in  the  region  just  above  the 
boundary  of  the  NR  region.  In  other  words,  it  was 
assumed  that  the  NR  boundary  line  corresponds  to 
the  line  “w  =  0°.  At  higher  Mach  numbers  Ms  4, 
the  line  0W  =  0°  (strictly,  the  line  which  is  the 
limit  of  *>w  -*•  0°)  coincides  with  the  NR  boundary 
line  within  the  error  of  calculations.  At  lower 
Mach  numbers,  however,  the  line  -  q°  separates 
from  the  boundary  of  the  NR  region.  Consequently, 
there  arises  a  region  between  the  ow  =  0°  line  and 
the  boundary  of  the  NR  region  in  which  no  three- 
shock-theory  solution  exists  (that  is,  no  physical¬ 


ly-realistic  solution  can  exist  even  if  solutions 
with  negative  0W  are  possible).  The  non-solution 
region  for  a  gas  with  constant  >  =  1.4  is  shown  in 
Fig.  A-3 . 

Experimentally,  many  points  exist  in  this  no¬ 
solution  region.  Consequently,  the  theory  is 
inadequate  in  this  region.  The  existence  of  the 
no-solution  region  corresponds  to  the  fact  that  > 
continues  to  increase  as  Ms  -  1 ,  as  described  in 
Appendix  B . 1 . 


Since  the  existence  of  this  no-solution  region 
is  associated  with  the  problem  in  Appendix  B.l,  the 
resolution  of  this  unrealistic  region  might  be 
possible  in  an  analysis  which  takes  into  account 
the  effect  of  the  bow  shock,  as  described  in 
Appendix  B.l. 

B . 3  Existence  of  Another  CMR  Region 

Another  strange  result  due  to  the  three-shock 
theory  is  that  it  indicates  the  existence  of  another 
CMR  region  near  the  NR  boundary.  In  this  section, 
the  usual  transition  criterion  between  SMR  and  CMR, 
Myr  <  1,  is  adopted.  When  viewing  the  variation 
of  M2T  with  "w  at  a  fixed  Ms  =  1.8  (Fig.  A-4 )  for  a 
gas  with  constant  ,  =  1.4,  for  example,  Myj  is  1.081 
at  w  =  4S°.  It  means  that  the  point  P  at  Ms =  1.8 
and  -w  =  45°  lies  in  a  CMR  region.  The  value  of 
Mjp  decreases  as  "u,  decreases.  Mvp  is  sonic  at 
"K  »  39.3 2°  and  then  enters  the  SMR  region.  The 
decreasing  value  of  Mjt  is  due  to  two  causes.  One 
is  that  Mj  |-  decreases  with  decreasing  C'w  and  the 
other  is  that  :j  approaches  90°  with  decreasing  w . 
The  angle  4j  becomes  90°  at  about  6W  =  13°  and 
keeps  increasing  with  decreasing  “w  to  values 
greater  than  90°.  An  increase  in  tj  90°  has  an 
effect  to  increase  M2T  contrary  to  1 ]  <  90°.  At  a 
certain  '-w,  this  effect  overrides  the  decreasing 
effect  due  to  the  decrease  of  Mjt,  and  M2T  berins 
to  increase.  Eventually,  M2T  becomes  sonic  again 
at  about  ow  =3°.  At  wedge  angles  below  this  value 
of  "w ,  M2T  is  greater  than  1,  which  means  that  it 
is  in  another  CMR  region  according  to  the  usual 
criterion . 

This  CMR  region  in  the  case  of  a  gas  with 
constant  ,  =  1 .4  is  shown  in  Fig.  A-S.  This  CMR 
region  is  not  realistic  because  the  validity  of 
the  solutions  of  the  three-shock  theory  is  in  doubt 
in  the  region  of  small  ;,w  and  low  Ms,  as  described 
in  Appendices  B.l  and  B.2.  This  additional  CMR 
region  does  not  appear  when  the  new  condition 
‘  ■  90°  (see  Section  3.3)  is  added  *o  the  criterion 
because  A  in  this  region  is  well  below  90°. 


APPENDIX  C 


EFFECT  OF  SLIPSTREAM  THICKNESS 


Only  the  flow  direction  and  pressure  are  iden¬ 
tical  in  regions  2  and  3  separated  by  the  slipstream 
but  not  the  other  physical  quantities.  Consequently, 
a  thermal  and  velocity  layer  results,  not  unlike  a 


shear  or  boundary  layer.  It  can  be  seen  from  Figs. 

5  and  A-10  (and  other  photographs  in  Refs.  1,  3,  4, 

13  and  14)  that  this  is  an  ideal  way  of  generating  a 
shear  layer  through  wave  interactions.  The  layer  is 


I 


APPENDIX  A 

COMMENT  ON  PERFECT  CO; 


The  term  perfect  CO2  was  used  previously  in  Refs. 
13  and  14  to  indicate  a  gas  with  a  constant  value  of 
V  =  1.29  independent  of  temperature.  However,  from 
the  following  discussion,  CO2  with  a  constant 
1  =  1.29  has  no  physical  validity.  At  room  tempera¬ 
ture,  the  vibrational  modes  of  CO2  are  partly  ex¬ 
cited.  The  doubly-degenerate  bending  modes,  which 
have  a  characteristic  temperature  of  960.2  K,  are 
excited  to  14%  of  their  full  excitation  and  have  a 
contribution  of  0.9R  to  Cv-  The  other  two  stretch¬ 
ing  modes  contribute  0.7R  to  Cv-  Consequently, 
along  with  contributions  from  the  translational 
modes  and  rotational  modes  Cv  =  3.47R  at  room  tem¬ 
perature,  or  1  =  Cp/Cv  =  1.29  at  room  temperature. 
However,  behind  a  shock  wave  the  translational  and 


rotational  temperatures  immediately  jump  to  the 
Rankine-Hugoniot  temperature  Tly.  Since  the  vib¬ 
rational  energy  is  too  slow  to  participate,  it 
remains  frozen  on  either  side  of  the  shock  front 
and  the  temperature  Tip  can  be  computed  using 
Y  =  1.4  (Ref.  31).  After  a  period  (relaxation 
time)  the  vibrational  modes  are  activated  and 
share  in  the  energy  through  equipartition.  The 
final  equilibrium  temperature  Tle  <  Tlf.  Perfect 
CO2  with  >  =  1.29  has  no  physical  basis.  However, 
it  is  sometimes  used  in  engineering  applications 
to  give  an  approximate  answer  for  shock  waves  in 
CO 2  with  vibrational  excitation.  The  error 
berrmes  increasingly  worse  with  rising  shock 
M  .lumber  (Ref.  31). 


APPENDIX  B 

THREE -SHOCK  THEORY 


Ihe  three-shock  theory  provides  a  simple  and 
effective  analysis  of  obi ique -shock -wave  reflec¬ 
tions.  Even  so,  it  includes  some  aspects  which 
are  physically  unreasonable  or  inconsistent  with 
experimental  results.  They  are  described  below. 

B.l  Behaviour  of  .  as  M  -*1 
- - -  s  - 

The  values  of  .  at  fixed  'w  are  approximately 
constant  at  high  shock  Mach  numbers  Ms,  both 
experimentally  and  analytically.  As  Ms  decreases, 
the  values  of  .  increase  except  for  large  9W. 
Experimentally,  .  has  maximum  values  near  Ms  =  2 
and  decreases  for  other  Ms.  The  calculated  results 
of  .  from  the  three-shock  theory,  however,  do  not 
have  maxima  for  different  but  continue  to  in¬ 
crease  az  M-  -l,  as  shown  in  Fig.  A-l. 

Ben-Dor  |Rrf.  3|  proposed  an  alternative  method 
to  obtain  .  at  small  wedge  angles  (9W  <  4°)  which 
does  give  maxima  and  may  agree  with  experiment 
despite  his  failure  to  obtain  solutions  of  the 
three-shock  ti.?ory  for  small  Hw  and  low  Ms .  His 
method  is  to  assume  that  at  small  0W  coincides 
with  the  value  at  the  no  reflection  (NR)  boundary 
’wNR  and  >  -  ’wNR  "  w  at  small  wedge  angles.  Des¬ 
pite  that  his  method  gives  a  reasonable  approxima¬ 
tion  to  experiment  results  and  is  useful  from  an 
engineering  point  of  view,  he  gave  no  physical 
explanation  that  at  small  coincides  with  the 
NR  boundary.  His  method  does  not  give  the  neces¬ 
sary  dependency  of  «w  (*  PwNR  >  for  different  Gw. 
Since  his  suggestion  is  outside  the  three-shock 
theory  realm,  being  only  an  assumption,  it  does 
not  provide  any  clue  for  improving  the  three-shock 
theory  or  an  alternative  theory  for  the  calculation 
of  / .  The  tendency  to  increase  .  as  Ms  *  1  is  not 
considered  realistic,  not  only  because  it  does  not 
agree  with  experimental  results,  but  it  also  yields 
a  no-solution  region  in  the  (Ms-0w)  plane,  as  des¬ 
cribed  in  Appendix  B.2. 


The  reason  why  the  three-shock  theory  fails  to 
give  reasonable  solutions  at  low  Mach  numbers  may 
be  as  follows.  The  reflected  shock  wave  is  com¬ 
posed  of  two  parts,  one  is  a  reflected  shock  wave 
emanating  from  the  triple  point  and  the  other  is 
a  bow  shock  wave  caused  by  the  interaction  of  the 
incident  flow  with  the  flow  deflected  by  the  wedge. 
The  two  shock  waves  intersect  at  the  kink  or  the 
second  triple  point  in  CMR  or  DMR ,  respectively. 

In  SMR,  the  two  shock  waves  blend  smoothly  so  that 
an  intersection  point  cannot  be  seen.  The  three- 
shock  theory  gives  the  position  of  only  one  of  the 
two  shock  waves,  which  is  the  one  emanating  from 
the  triple  point.  The  actual  shape  of  the  entire 
reflected  shock  wave  is  considered  to  be  determined 
by  the  two  shock  waves  mentioned  above  and  their 
interaction.  At  low  Mach  numbers,  both  the  reflect 
ed  shock  wave  from  the  triple  point  and  the  bow 
shock  wave  are  weak.  The  former,  however,  is  con¬ 
sidered  to  become  weaker  more  rapidly  than  the 
latter.  Therefore,  the  position  of  the  whole  re¬ 
flected  shock  wave  is  controlled  by  the  bow  shock 
wave  rather  than  the  reflected  shock  wave  from 
the  triple  point.  In  the  limit  as  the  shock  Mach 
number  of  the  bow  shock  wave  approaches  1  (Mach 
wave),  or  -  0°,  the  bow  shock  will  intersect 
the  incident  shock  wave  at  a  point  =  wnr.  which 
is  at  a  lower  0^  than  the  limit  of  solution  of  the 
three-shock  theory. 

The  reason  is  illustrated  in  Fig.  A-2.  Point 
0  in  Fig.  A-2  was  at  the  wedge  comer  when  the 
incident  shock  wave  collided  with  the  wedge  corner 
and  it  moves  with  the  flow  in  region  1_.  Therefore, 
the  flow  in  region  1  has  a  direction  TO  relative 
to  T  and  the  bow  shock  wave  in  the  vicinity  of  T  is 
normal  to  the  flow.  Since  the  bow  wave  is  a  Mach 
wave,  the  flow  velocity  in  region  1  relative  to  the 
bow  wave,  or  T,  is  sonic.  If  T  is  at  ft^NR,  the 
flow  Mach  number  in  region  1  is  also  1.0  according 
to  the  criterion  of  the  NR  boundary.  It  means  that 


/ 


laminar  near  its  point  of  generation,  the  triple 
point,  and  further  on  it  becomes  turbulent  and 
thicker.  The  thickness  of  the  slipstream  between 
the  two  flows  violates  one  of  the  boundary  condi¬ 
tions  of  the  three-shock  theory  that  the  two  flows 
are  parallel.  In  this  appendix,  the  actual  dis¬ 
placement  thickness  of  the  slipstream  between 
regions  2  and  3  is  not  estimated.  Only  the  effect 
of  an  assumed  displacement  angle,  which  is  the 
slope  of  the  displacement  thickness  at  the  triple 
point  on  \  is  estimated  in  order  to  evaluate  whether 
or  not  this  effect  can  explain  the  remaining  dis¬ 
crepancies  of  x  between  the  experiments  and  the 
analytical  predictions. 

Figure  A-6  shows  a  comparison  of  calculated 
results  with  and  without  displacement  angles  of  1° 
and  2°.  The  gas  is  equilibrium  CO2  at  Ms  =  5. 


Since  the  figure  shows  the  relation  between  bw  and 
9*.  the  difference  between  the  value  along  the 
longitudinal  and  transverse  axis  is  the  value  of  x- 
A  displacement  angle  of  2°  can  change  the  value  of 
X  less  than  1°,  except  at  small  wedge  angles 
(6w  <  10°).  This  is  insufficient  to  explain  the 
discrepancy  in  x  (which  is  about  2°  in  a  typical 
case  at  Ms  =  6  and  9W  =  20°)  because  a  displacement 
angle  of  6°  is  necessary  to  explain  the  discrepancy 
and  it  is  not  likely  to  happen.  It  should  be 
pointed  out  here,  however,  that  the  effect  of  the 
displacement  angle  is  not  negligibly  small  at 
small  wedge  angles.  The  displacement  angle  should 
be  taken  into  account  in  calculations  at  small 
wedge  angles  regardless  whether  use  is  made  of  the 
three-shock  theory  which  takes  into  account  the 
effect  of  the  bow  shock  wave,  as  suggested  in 
Appendix  B.l. 


APPENDIX  D 


COMMENT  ON  TRANSITION  CRITERION  FOR  CMR  TO  DMR 


The  discrepancy  between  the  calculated  CMR  *  DMR 
boundary  and  experiment  is  not  as  simple  as  for  Che 
SMR  t  CMR  boundary.  In  the  case  of  COz  shown  in 
Fig.  18,  the  calculated  boundary  fails  to  account 
for  many  CMR  points  which  lie  in  the  DMR  region. 

That  is,  the  CMR  *  DMR  boundary  line  must  be  shifted 
upward  to  account  for  the  experimental  results, 
whereas  in  the  case  of  Ar  shown  in  Fig.  A-7,  the 
reverse  occurs  and  the  CMR  *  DMR  line  has  to  be 
shifted  downward.  However,  for  Nj  and  air,  which 
are  shown  in  Figs.  A-8  and  A-9,  respectively,  the 
experimental  distributions  agree  well  with  the 
calculated  boundaries  (except  at  lower  Ms  and  large 
9W  when  the  two  triple  points  approach  and  coincide 
on  the  RR  boundary  line).  All  calculations  were 
done  in  accordance  with  the  discussion  on  relaxation 
lengths,  that  is,  C02  us  an  equilibrium  gas  and  Ar, 
Nz  and  air  as  a  frozen  gas.  These  discrepancies  are 
considered  due  to  at  least  two  causes.  One  is 
related  to  the  classification  of  the  experimental 
results  and  the  other  is  due  to  the  inappropriate 
transition  criterion. 

Experimentally,  the  distinction  between  CMR 
and  DMR  is  made  by  whether  or  not  there  is  another 
shock  wave  at  the  second  triple  point,  it  is 
difficult,  however,  to  see  a  weak  and  short  shock 
wave  in  an  interferogram.  Consequently,  the  classi¬ 
fication  sometimes  depends  on  the  observers.  For 
example,  the  interferogram  of  Ar  shown  in  Fig.  A-10 
was  classified  as  a  DMR  by  Ben-Dor  (Ref.  4)  and  the 
one  of  COj  shown  in  Fig.  A-ll  was  classified  as  a 
CMR  by  Ando  (Ref.  14).  There  seems,  however,  no 
essential  difference  between  these  two  interfero- 
grams,  or  one  can,  with  greater  assurance,  classify 
them  just  in  reverse,  that  is.  Fig.  A-10  as  a  CMR 
and  Fig.  A-ll  as  a  WR.  If  so,  the  experimental 
boundary  of  Ar  may  jove  towards  the  DMR  region  and 


that  of  CO2  towards  the  CMR  region. 

From  the  point  of  the  criterion,  a  criticism 
described  below  can  be  offered.  At  a  flow  Mach 
number  just  over  unity,  a  shock  wave,  if  it  is 
formed,  must  be  normal  to  the  flow  direction,  other¬ 
wise  the  Mach  number  of  the  normal  component  of  the 
flow  is  less  than  1  and  no  shock  wave  could  be 
formed.  As  mentioned  in  Section  3.3,  the  flow  comes 
right  above  (parallel  to  the  incident  shock -wave 
front),  looked  at  from  the  reference  frame  attached 
to  the  second  triple  point.  In  the  region  where 
the  transition  from  CMR  to  DMR  takes  place,  the 
angle  i  is  larger  than  90°  (see  Fig.  20).  Therefore, 
the  flow  will  be  deflected  toward  the  second  triple 
point  T'  while  passing  through  the  reflected  shock 
wave.  As  a  result,  the  second  Mach  stem  M'  must 
slant  below  the  horizontal  line  toward  the  Mach  stem 
M  (see  Fig.  A- 12).  The  shadowgraph  shown  in  Fig.  17 
is  just  inside  the  DMR  region.  The  direction  of  the 
flow  after  the  deflection  through  the  reflected 
shock  wave,  looked  at  from  the  second  triple  point, 
is  calculated  at  an  angle  of  about  5°  toward  the 
second  triple  point  from  the  vertical  line.  There¬ 
fore,  the  second  Mach  stem  should  lie  5°  below  the 
horizontal  line  toward  the  Mach  stem.  Experimentally, 
however,  the  second  Mach  stem  lies  along  the  hori¬ 
zontal  line,  or  it  seems  to  have  an  opposite  inclin¬ 
ation  in  the  vicinity  of  the  second  triple  point. 

The  second  Mach  stem  seems  to  be  formed  as  an  exten¬ 
sion  of  the  bow  shock  wave  rather  than  as  a  shock 
wave  normal  to  the  flow  in  region  2  at  the  very 
beginning  of  its  formation.  This  difference  is, 
however,  too  small  to  explain  the  whole  discrepancy 
between  the  experiments  and  the  calculated  boundary, 
especially  in  the  case  of  CO2 .  It  gives  only  a  clue 
to  other  necessary  conditions  for  the  transition 
lines  which  may  exist. 


APPENDIX  E 


SOME  FUTURE  STUDIES 


There  are  still  many  aspects  to  be  studied  ex¬ 
perimentally  and  analytically  in  the  area  of  oblique- 
shock-wave  reflections.  Some  of  these  are  summarized 
in  this  appendix  in  order  to  assist  in  the  continua¬ 
tion  of  the  study  of  oblique-shock-wave  reflections. 

(a)  Experimental  persistence  of  RR  beyond  the  tran¬ 
sition  boundary. 

More  experiments  are  required  prior  to  further 
analyses.  Although  the  persistence  of  RR  beyond  the 
transition  boundary  of  the  detachment  criterion 
itself  is  an  apparent  experimental  fact,  some  ex¬ 
perimental  results  seem  inconsistent  with  each 
other.  That  is,  both  RR  and  MR  are  obtained  in  the 
same  region  of  the  tMs-Ow)  plane  from  the  results 
of  different  experiments  (see  Section  3.4).  This 
might  be  due  to  unknown  factors  such  as  a  difference 
in  initial  pressure  and  its  effects  on  the  boundary- 
layer-displacement  angle  Dj.  A  knowledge  of  the 
exact  position  of  the  experimental  transition  boun¬ 
dary  and  its  dependence  on  controlled  initial 
conditions  are  required  for  further  analyses  of  this 
problem.  The  experimental  data  are  required  not 
only  at  the  boundary  but  also  above  and  below  the 
boundary  to  obtain  information  on  the  behaviour  of 
quantities  such  as  '  (in  RR)  and  x  (in  MR)  when 
approaching  the  RR  t  MR  boundary  line. 

(b)  Convergence  of  the  SMR  t  CMR  and  CMR  *  DMR 

boundaries  at  the  RR  *  MR  boundary. 

Experimental  results  show  that  the  CMR  *  DMR 
boundary  approaches  the  SMR  *  CMR  boundary  near 
the  RR  boundary  and  they  eventually  converge  at 
the  RR  *  MR  boundary  (see  Fig.  25).  This  corre¬ 
sponds  to  the  experimental  fact  that  V  L  approaches 
0  and  0w  comes  near  to  the  RR  *  MR  boundary.  There 
is  no  analytical  prediction  of  the  position  of  the 
second  triple  point  except  for  an  empirical  assump¬ 
tion  that  t/L  is  equal  to  the  ratio  between  the 
velocity  of  the  flow  behind  the  incident  shock  wave 
and  the  incident  shock  wave  velocity  ui^/usn.  An 
analytical  method  to  predict  the  position  of  the 
second  triple  point,  especially  near  the  RR  boun¬ 
dary,  is  required.  Experiments  are  also  required 
since  few  data  are  available  aside  from  the  existing 
values  of  i/l  in  this  region. 


(c)  Prediction  of  the  direction  of  the  Mach-stem 

curvature . 

As  described  in  Section  3.3,  the  Mach  stem  at 
the  tiple  point  is  not  perpendicular  to  the  wedge 
surface  experimentally.  This  is  because  the  Mach 
stem  has  curvature.  This  discrepancy  affects  the 
analytical  solutions,  especially  the  value  of  y . 

For  accurate  analyses,  the  direction  of  the  Mach 
stem  at  the  triple  point  must  be  predicted.  A 
semi-empirical  formulation  will  be  helpful  until 
a  better  solution  is  obtained,  which  requires 
solving  the  two-dimensional  flow  field  itself 
without  numerically  smearing  the  shock  fronts. 

(d)  Effect  of  the  bow  shock  wave  on  the  position 

of  the  triple  point. 

As  described  in  Appendix  B.l,  the  three-shock 
theory  fails  to  predict  x  at  low  Ms  and  small  6W. 

This  is  believed  to  be  due  to  the  effect  of  the  bow 
shock  wave  (see  Appendix  B.l).  An  alaytical  method 
to  predict  the  strength  of  the  bow  shock  wave  and 
its  effect  on  the  position  of  the  triple  point  is 
required . 

(e)  Transition  boundaries  in  the  MR  region. 

There  is  experimental  impreciseness  in  the  tran¬ 
sition  boundaries  for  the  types  of  MR.  The  impreci¬ 
sion  in  the  classification  between  SMR  and  CMR 
arises  from  the  definition  of  CMR  itself.  In  the 
case  between  CMR  and  DMR,  however,  an  interferogram 
(which  responds  to  changes  in  r.  unlike  schlieren 
to  di./dx  and  shadowgraph  to  d2c/dx2)  contributes  to 
the  uncertainty  of  the  experimental  transition  line. 

A  shock  wave  is  easier  to  distinguish  in  a  shadow  or 
schlieren  photograph  rather  than  in  an  interferogram 
(as  noted  above) .  It  is  true  that  an  interferogram 
includes  more  quantitative  information  than  a 
shadowgraph.  Even  so,  experiments  using  shadowgraphs 
are  often  required  to  give  precise  configurations 
which  determine  the  positions  of  the  experimental 
transition  boundaries,  especially  between  CMR  and 
DMR.  Consequently,  more  than  one  optical  method  of 
flow  visualization  should  be  used. 


FIG.  A-l  TRIPLE-POINT-TRAJECTORY  ANGLE  x  VS  INCIDENT  SHOCK  MACH  NUMBER  Ms 
FROZEN  GAS,  y  =  1.40. 
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FIG.  A-5  EXISTENCE  OF  ANOTHER  CMR 
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FIG.  A-8  COMPARISON  OF  CALCULATED  CMR  t  DMR  TRANSITION  BOUNDARY  WITH 
EXPERIMENTAL  RESULTS  FOR  N2 . 

EXPERIMENTS:  o  -  RR,  A  -  SMR,  V  -  CMR,  Q  -  DMR  (REF.  4). 
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FIG.  A-9  COMPARISON  OF  CALCULATED  CMR  t  DMR  TRANSITION  BOUNDARY  WITH 
EXPERIMENTAL  RESULTS  FOR* AIR. 

EXPERIMENTS:  o  -  RR,  A  -  SMR,  V  -  CMR,  o  -  DMR  (REF. 21  ). 
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FIG.  A-12  SCHEMATIC  DIAGRAM  ILLUSTRATING  THE  DIRECTION  OF  SECOND  MACH  STEM 
AT  ITS  VERY  FIRST  APPEARANCE. 

-  SHOCK  WAVE, - SLIPSTREAM,  — *■  FLOW  DIRECTION 

(RELATIVE  TO  T’),  -  HORIZONTAL  LINE,  I  -  INCIDENT  SHOCK 

WAVE,  R,  R'  -  REFLECTED  SHOCK  WAVE,  M  -  MACH  STEM,  M’  -  SECOND 
MACH  STEM,  T  -  TRIPLE  POINT,  T'  -  SECOND  TRIPLE  POINT,  S  -  SLIP- 
STREA,  S'  -  SECOND  SLIPSTREAM. 
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